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Abstract. In this paper the spin configurations of the ground state and one- and 
two-electron excited states of the Hubbard model on the square lattice are studied. 
^ ' We profit from a general rotated-electron description, which is consistent with the 

model global SO{3) x SO{3) x C/(l) symmetry. For rotated electrons, doubly and 
^ ' single occupancy are good quantum numbers for on-site repulsion U > 0. The above 

states are within that description generated by occupancy configurations of charge c 
fermions and spin-singlet two-spinon si bond particles. Those describe the charge and 
spin degrees of freedom, respectively, of the rotated electrons that singly occupy sites, 
i-rt ■ While the c fermions have no internal structure, that of the spin-neutral si bond- 

^ ' particle occupancy configurations is here described in terms of spinon occupancies of 

a well-defined effective spin lattice. In reference [2] it is confirmed that our results 
contribute to the further understanding of the role of electronic correlations in the 
^^ spin spectrum of the parent compound La2Cu04. They are also of interest for studies 

^ . of ultra-cold fermionic atoms on an optical lattice. 
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1. Introduction 

The Hubbard model on a square lattice is the simplest toy model for the description 
of the effects of electronic correlations in the spin spectrum of the Mott-Hubbard 
insulators cuprates parent compounds [H |2]. On the square and cubic lattices the 
model has no exact solution and many open questions about its properties remain 
unsolved. Fortunately, it can be experimentally realized with unprecedented precision 
in systems of correlated ultra-cold fermionic atoms on an optical lattice [3]. Recent 
experiments considered the Hubbard model on the cubic lattice [1]. They involved 
ultra-cold fermionic atoms on an optical cubic lattice formed by interfering laser fields. 
Similar experiments referring to the model on the square lattice are in progress. One 
may expect very detailed experimental results over a wide range of parameters to be 
available. 

In this paper we study the spin occupancy configurations of the ground state of 
the Hubbard model on a square lattice with A^^ ^ 1 sites within the spin-singlet two- 
spinon si bond-particle description also used in the investigations of Ref. [2]. Here Na 
is the number of sites of an edge of length L = Naa and spacing a. Hence the total 
number of sites of the square lattice of area L^ is A^^ = Na x Na. That description 
involves the generalization to all finite values of the on-site repulsion f/ > of the exact 
transformation for separation of spin-1/2 fermions without constraints introduced for 
large U values in Ref. [5]. It is consistent with the global S0{3) x 5*0(3) x f/(l) = 
[5*0(4) X U{l)]/Z2 symmetry found recently in Ref. p] for the model on any bipartite 
lattice. That global symmetry is an extension of the 50(4) symmetry known to occur for 
the model on such lattices |7]. The extended global symmetry is related to the rotated 
electrons that for U > emerge from the electrons through a unitary transformation of 
the type considered in Ref. [8] and to the local symmetries and unitary transformations 
considered in Ref. [9]. For rotated electrons double and single occupancy are good 
quantum numbers for U > 0. 

Here we study the spin configurations of the model ground states and excited states 
that span the subspace relevant for the one- and two-electron physics. It is spanned by 
an initial vanishing-spin-density ground state and the energy and momentum eigenstates 
contained in the excitations generated by application onto that ground state of one- and 
two-electron operators. Within the description used in the studies of Ref. [2], for N^ ^ 1 
such states can be generated by suitable occupancy configurations of two basic quantum 
objects: spin-less and ?7-spin-less charge c fermions and spin-1/2 spinous. There are also 
?7-spin-l/2 //-spinous, yet in that subspace they refer to a single occupancy configuration 
and thus play no active role. The spinous describe the spin degrees of freedom of the 
rotated-electron singly occupied sites. Their spin projection is the same as that of the 
corresponding rotated electron. The ?7-spinons of ?7-spin projection down and up describe 
the 77-spin degrees of freedom of the rotated-electron doubly occupied and unoccupied 
sites, respectively. In our subspace rotated-electron double occupancy vanishes and 
there are only ?7-spin-up //-spinous. The spinon spin-singlet configurations involve in 
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general the 2z/-spinon and 2z/-site su bond particles [2]. Here z/ = 1, 2, 3, ... refers to the 
number of spinon pairs in each su bond. Fortunately, the ground states and their one- 
and two-electron excited states considered here, except the spin-singlet states, involve 
only two-spinon and two-site si bond particles. Moreover, the latter states involve in 
addition to the si bond particles a single four-spinon and four-site s2 bond particle. 

The c fermions live on a c effective lattice identical to the original lattice. The c 
fermion occupied sites describe the charge degrees of freedom to those singly occupied 
by rotated electrons. The corresponding spin degrees of freedom are described by the 
spinous. In turn, the c fermion unoccupied sites correspond to those doubly occupied 
or unoccupied by rotated electrons. The corresponding ry-spin degrees of freedom are 
described by the //-spinous. Importantly, the spatial coordinates of the c fermion 
occupied and unoccupied sites are also the spatial coordinates in the original lattice 
of the sites referring to spinous and ?7-spinons, respectively. Provided that N^ ^ 1 this 
is behind the occupancies of the spin effective lattice and ?7-spin effective lattice being 
independent [2]. Such effective lattices involve only the rotated-electron singly occupied 
sites and the rotated-electron doubly occupied and unoccupied sites, respectively. 

The spinous that are not part of su bond particles are called independent spinous. 
They are invariant under the electron - rotated-electron unitary transformation. For 
the lowest-weight states whose spin Sg and spin projection S^ are such that Sg = —S^ 
the independent spinous have up spin projection. For the subspace considered here all 
?7-spinons have up 77-spin projection so that the r^-spinon effective lattice plays no active 
role. In turn, since the location in the original lattice of the sites of the spin effective 
lattice is recorded in the c fermion occupancy configurations, it turns out that provided 
that N^ :^ 1 it is a good approximation to consider that the spin effective lattice is a 
square lattice with the same length edge L as the original lattice and a number of sites 
A^^^ and spacing a^ given by, 

N^^ = il-x)N^; a, = ^^-^, (1 - x) > l/N^ . (1) 

Here D = 2 refers to the present square lattice, D = 1 corresponds to the one- 
dimensional (ID) lattice some times considered below, and x = [N^ — N)/N^ is 
the hole concentration. For A*"^ ^ 1, in spite of rotated-electron single occupancy 
being only conserved globally and not locally, the spinon occupancies of the spin 
effective lattice provide a good approximation of the ground state and excited states 
spin configurations [2]. Furthermore, within the present description such configurations 
involve the occupancies of a si effective lattice. Its unoccupied sites are the up-spin 
independent spinous. Its occupied sites refer to the two-spinon si bond particles studied 
in this paper. Hence each occupied site of the si effective lattice corresponds to two 
sites of both the original lattice and spin effective lattice. 

Bond states based on two-spinon and two-site bonds are in general not orthogonal 
and their basis is overcomplete. However, within the representation used in the studies 
of this paper such a problem does not occur. This is due to physical restrictions of 
the occupancy configurations of the well-defined set of two-site bonds contributing to 
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a given si bond particle. For instance, that they are centered at the same real-space 
coordinate of the spin effective lattice. For ID the discrete momentum values of both the 
c fermions and si fermions generated from the si bond particles studied here through 
an extended Jordan- Wigner transformation [21 [TOl [TTl [12] are good quantum numbers 
whose occupancy configurations generate the energy eigenstates [13]. 

Evidence that for the one- and two-electron subspace the states generated by the 
occupancy configurations of such discrete momentum values are for the model on the 
square lattice energy eigenstates is given in the related investigations carried out in 
Ref. 0. The square-lattice quantum liquid introduced in that reference contains the 
one- and two-electron excitations of the Hubbard model on a square lattice. At hole 
concentration x = {N^ — N)/Nl = 0, U/At ^ 1.525, and t ^ 295 meV it is found 
in that reference to quantitatively describing the spin-wave spectrum observed in the 
parent compound La2Cu04 [Ij. In turn, it is expected that the description of the role 
of electronic correlations in the unusual properties of the cuprate hight-temperature 
superconductors [TH [151 [iHl [E] requires as well accounting for the effects of three- 
dimensional uniaxial anisotropy and intrinsic disorder. 

The spin-singlet two-spinon si bond particle configurations studied in this paper are 
an example of the general resonating-valence-bond pictures for spin-singlet occupancy 
configurations of ground states studied in Refs. [HI [19]. As mentioned above, the 
particular type of spin configurations based on two-spinon and two-site bonds considered 
in our studies lack both the non-orthogonality and overcompleteness problems. Such 
spin-singlet two-spinon si bond particle configurations refer to the c and si fermion 
description also used in the investigations of Ref. [2] on the Hubbard model on the 
square lattice. Some progress in understanding the physics of that model has been 
achieved for different limits, through a variety of methods. It corresponds to a non- 
perturbative quantum problem in terms of electron operators, so that rewriting the 
theory in terms of the standard formalism of many-electron physics is an extremely 
complex problem. A detailed and extensive discussion of the relation between the new 
results obtained both in Ref. [2] and this paper by means of the c and si fermion 
description and previously known results is presented in that reference. 

For instance, it is shown that the predictions of the square-lattice quantum liquid 
theory concerning the spin spectrum at half filling as described in terms of the spin- 
singlet two-spinon configurations studied in this paper agree both with experiments 
on the parent compound La2Cu04 and results obtained by the standard formalism of 
many-body physics. Furthermore, in Ref. [2] the relation of the description used in the 
studies of this paper and that reference to several other schemes is discussed. 

The paper is organized as follows. The model, the basic rotated-electron 
representation, and the description in terms of the basic quantum objects, which emerge 
from the rotated-electron occupancy configurations, are the subjects of Section 2. The 
si effective lattice, corresponding change of gauge structure, and si operators of the 
N^i = configuration state are the topics addressed in Section 3. (Here A^^^^ denotes 
the number of unoccupied sites of the si effective lattice.) Section 4 contains a study 
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of the N^^ = 1,2 configuration states and kink-like and anti- kink-like link occupancies 
associated with the si effective lattice unoccupied sites. Finally, Section 5 contains the 
concluding remarks. 

2. The Hubbard model on the square lattice, rotated electrons, and the 
description used in our studies 

2.1. The Hubbard model on the square lattice and rotated electrons 

The Hubbard model on the two-dimensional (2D) square lattice with torus periodic 
boundary conditions and the same model on the ID lattice with periodic boundary 
conditions, spacing a, N^ = [Na]^ sites where D = 1 and Z^ = 2 for the ID and square 
lattices, respectively, N^ ^ 1 even, and lattice edge length L = Naa for 2D and chain 
length L = A'^a a for ID is given by, 

H = -tY. E l^l^ ^'--^ + ^•^•] + U i^^ - ^]/2 • (2) 

Here the operator, 

Q = ^^ ^r„a (1 - n^„-a) , (3) 

where n^-. o- = ct ^c^^^a and —a =t (and — cr =1) for a =J, (and a =t) counts the number 
of electron singly occupied sites. Hence the operator D = [N — Q]/2 counts the number 
of electron doubly occupied sites where N = ^^ N^ and iVg- = J2j=i ^rj,a- 

The studies of Ref. [2] use a uniquely-defined electron - rotated-electron unitary 
V = V{U/At). It is such that states |^t//4i) = V^'^l^oo) are for U/At > energy 
eigenstates. It corresponds to a suitable chosen set {|\l/oo)} of U/At — )■ oo energy 
eigenstates. The unitary transformation maps the electronic operators ct and c^ ^ 
onto rotated-electron creation and annihilation operators ct = yt 4 y and c^. „ = 

^ rj,a rj,a 'jt" 

V^ Cf^j,a V, respectively. 

The studies of Ref. [6] reveal that for U/At > the Hubbard model on a square 
lattice has a global 5*0(3) x SO (3) x f/(l) symmetry. The generator Sc of the hidden 
global U{1) symmetry reads Sc = V'^ ScV where Sc = Q/2 and the operator Q is given 
in Eq. ([3]). Its eigenvalue Sc is one-half the number of rotated-electron singly occupied 
sites 2Sc. Indeed, the rotated electrons of Refs. [21 [6] are constructed in such a way that 
Sc is a good quantum number in the rotated-electron picture. The subspaces spanned 
by states whose number 2Sc is constant play an important role. For hole concentrations 
< a; < 1 and maximum spin density m = (1 — a;) there is a fully polarized vacuum, 
which remains invariant under the electron - rotated-electron unitary transformation, 

10,3) = |0,; <) X |0,; O X \GSc; 2Sc) . (4) 

Here the r^-spin SU{2) vacuum |0^; Nj^) associated with A^^ = [A^^ — 2S'c] independent 
-M/2 ?7-spinons, the spin SU{2) vacuum |Os;A^jf) with A^^ = 2Sc independent -1-1/2 
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spinons, and the c f/(l) vacuum \GSc', 2Sc) with A^^ = 2S'c c fermions remain invariant 
under the electron - rotated-electron unitary transformation. In turn, for states with 
a finite number of si bond particles the c fermions are not invariant under that 
transformation. 

The whole physics can be extracted from the model ([2]) in the LWS-subspace 
referring to values of Sa and 5*^ such that Sa = —S^ for a = 77, s [2]. Within the 
LWS representation, the c fermion creation operator can be expressed in terms of the 
rotated-electron operators as follows, 

where n^.^„ = c\..^c^-o- and e*'^'^'^ is ±1 depending on which sublattice site Tj is on. The 
three spinon local operators si and three r^-spinon local operators pi such that / = ±, 2 
and sj = s% ±isl. and »i = p% ±ip^ , respectively, where the Cartesian coordinates 

rj Tj Tj i Tj 1- J j '■ 'rj' 

X, y, z are often denoted in this paper by X\,X2-, ^3, respectively, are given by, 

4- = ^r,,c (if, ; pt- = (1 - ^r, ,c) g'- , ^ = ±, ^ ; n^- ,c = f\^c fr,,c ■ (6) 

Here n^^^c = flcf^ ,c is the c fermion local density operator and the rotated quasi- 
spin operators read as follows in terms of rotated-electron creation and annihilation 
operators, 

gt = (gj.^^ - e'"'-''^ c,^^,t) cf,,i ; qz^ = (gt )^ ; qp, = i^- nf,A ■ (7) 



'J 



Inversion of the relations provided in Eqs. ([5]) and ([7]) gives, 

4,,, = %(4,c-e^""^4,c)- (8) 

For the one- and two-electron subspace N^ ^ 1 and the number A^^^ of sites of 
the si effective lattice, A^'^i of si bond particles, and N^i of unoccupied sites read, 

iViJ^ = iV,i + N^, = nZ/2 + Ss ; iV,i = nZ/2 - [S, + 2N,^] , 

N^, =[2S, + 2iV,2] = 0,l,2, (9) 

whereas for the corresponding c fermions such numbers are given by, 

NZ = N, + N^.=NJ^: N, = 2S, = {l-x)N^- N^ = x N^ . (10) 

For that subspace the si effective lattice is either full with A^^i = N^ = N^J2 or has 
one or two unoccupied sites. For it there is commensurability between the real-space 
distributions of the iV^^ ^ A^^i sites of the si effective lattice and N^^ ^ 2A'^si sites of the 
spin effective lattice. For (1 — x) > l/N^ the spin effective lattice has A*"^ = (1 — x) A^^ 
sites. The A''^^ expression given in Eq. Q) implies that a^i = L/Na^^ reads. 
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where 25"^ = 0,1,2 and the lattice constant a^ of the spin effective lattice is given in Eq. 

©• 

For A^j\ = states the bipartite ID and square spin effective lattices have two 
well-defined sub-lattices. For the square lattice the two spin effective sub-lattices have 
lattice constant a^i = v2as. In turn, for ID the sites of each spin effective sub-lattice 
are distributed alternately along the chain, the corresponding nearest-neighboring sites 
being separated by a^i = 2as. The fundamental translation vectors of such sub- lattices 
read, 

4i = ctsi el'i , [1^] ; flsi = -^(e^i-Fe^a) ; ^^i = — ^(ex^-ex^) , [2/^] , (12) 

where e^.^ and 6^2 ^^^^ the unit vectors. The vectors given in Eq. flT^ are the fundamental 
translation vectors of the si effective lattice [2]. 

The c fermions are ?7-spinless and spinless fermions without internal structure and 
their effective lattice is identical to the original lattice. In contrast, the composite two- 
spinon si bond particles have internal structure and expression of their si effective 
lattice occupancies in terms of spinon occupancies of the spin effective lattice and 
rotated-electron occupancies of the original lattice is a more complex problem, which 
deserves and requires further studies. This is the main subject of the remaining of this 
paper. 

3. The N!}-^ = configuration state: si effective lattice, corresponding 
change of gauge structure, and si bond-particle operators 

For the square- lattice model in the one- and two-electron subspace the numbers A^^i, 
A^^^, and A^g'\ = [A^i^^ — Ngi] are good quantum numbers whose values are fully 
controlled by those of S^, Sg, and Sc- We call configuration states the spinon occupancy 
configurations in the spin effective lattice which generate the spin degrees of freedom of 
the energy eigenstates that span the one- and two-electron subspace. Such configuration 
states refer to the overall occupancy of the A^^i si bond particles over the A*"^^ sites of 
the si effective lattice. Each configuration state refers to well-defined positions of the 
iV^i unoccupied sites. A x > and m = ground state is an example of an energy 
eigenstate whose spin degrees of freedom are described by the A^^^ = configuration 
state studied in the following. 

We start by considering the A^^^ = [A^^ — A^^i] = configuration state with no 
unoccupied sites in the spin and si effective lattices. It describes the spin degrees of 
freedom of x > and m = ground states and its charge excited states belonging to 
the one- and two-electron subspace whose A"^^ = N^J2 = (1 — x) N^/2 sites of the si 
effective lattice are occupied. Such a configuration state is described by suitable spinon 
occupancy configurations of the A"^ = (1 — x) A"^ sites of the spin effective lattice. 
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3.1. Independent two-site one-link bonds 

For simplicity we consider that Na^ = (1—xY'^ Na is an integer number yet if otherwise 
one reaches the same results in the N^ ^ 1 limit that our study refers to. For the model 
on the square lattice the spin effective lattice is then a square lattice with Na^ x A^^^ 
sites. For the model on that lattice we consider torus periodic boundary conditions for 
the spin effective lattice, alike for the original lattice. That implies periodic boundary 
conditions for the Na^ rows and A^^^ columns. Periodic boundary conditions are also 
used for the one-chain spin effective lattice of the ID model. 

The spin effective lattice has in the present case two sub-lattices. The real-space 
coordinates of the sites of each of such sub-lattices correspond to a possible choice of 
those of the si effective lattice. Indeed, there is for the N^-^ = configuration state a 
gauge "symmetry" between the representations in terms of the occupancies of the two 
alternative choices of real-space coordinates of the si effective lattice, which are found 
to refer to two alternative and equivalent representations of that state. We use the same 
notation as Xiao-Gang Wen in Ref. [22] and say that there is a gauge structure when 
we use simultaneously the real-space coordinates of the two corresponding choices of si 
effective lattices to label the N^^ = configuration state. The real-space coordinates 
fj of such two sub-lattices have j = 1,...,N^ sites, in either case the fundamental 
translation vectors being those given in Eq. (iT2l) . The real-space coordinates fj of the 
si bond particles are chosen to correspond to those of one of these two sub-lattices. 
Throughout the remaining of this paper we call sub-lattice 1 and sub-lattice 2 the sub- 
lattice of the spin effective lattice of a A^^'\ = configuration state whose real-space 
coordinates are and are not the same as those of the si effective lattice, respectively. 

An one-link bond connects two sites of the spin effective lattice whose spinous 
have opposite spin projection and correspond to a two-site spin-singlet configuration 
defined below. Each si bond particle is a suitable superposition of a well-defined set of 
two-site one-link bonds. Such a spin-singlet two-spinon si bond particle is related to 
the resonating- valence-bond pictures for spin-singlet occupancy configurations of ground 
states studied in Ref. [18], [19]. However, the si bond particle is well defined for all values 
of U /At > and not only for U/4t ^ 1, consistently with its two spinous referring to 
spins of the sites singly occupied by rotated electrons. In contrast, most schemes used 
previously for the Hubbard model or related models involving singly-occupied-site spins 
refer in general to large values of U/4t ^ 1 only [T6l [T2l [TSl [22] . Here the si bond 
particles have been constructed to inherently involving spinon occupancy configurations 
of sites of the spin effective lattice, which for U/At > refers only to the sites of the 
original lattice singly occupied by rotated electrons. 

For simplicity we often call two-site link or just link a two-site one-link bond. For 
the N^^ = configuration state studied here the above superposition includes 2D = 2,4 
families of two-site links, each family having Nsi/2D different types of such links: 
Nsi/2D is the largest number of independent links with the same link centre that exist 
for the above-considered boundary conditions. (Above and in the remaining of this 
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paper we denote often the number of family links by 2D = 2,4 where two and four is 
the number of such families for the model on the ID and square lattice, respectively.) 
Link independence means here that for a given link centre all links involve different pairs 
of sites and each site belongs to one pair only. 

The set of independent links with the same link centre belong to the same link 
family. Each one-link bond of a given family has some weight which in some cases may 
vanish. For a si bond particle of real-space coordinate r^ there are 2D = 2,4 families 
of links. The links of each family are centered at one of the 2D = 2,4 points fj + rj"; 
where the indexes d = 1,2 ioi D = 2, d = 1 for D = 1, and / = ±1 uniquely define 
the link family. Here f^i is the primary link vector. It connects the site of real-space 
coordinate fj to the centre of the four (and two for ID) links of real-space coordinate 
fj + f^i in the spin effective lattice involving that site and its nearest-neighboring sites. 
Note that the former site and the latter four (and two for ID) sites belong to sub-lattice 
1 and 2, respectively. On choosing one of the two sub-lattices of the spin effective lattice 
to be sub-lattice 1 and thus playing the role of si effective lattice and representing the 
states in terms of the occupancy configurations of the latter lattice we say that there is 
a change of gauge structure [22]. By considering all Nsi/2D independent links for each 
of the 2D link centres needed to describe a si bond particle of real-space coordinate fj 
we consider the most general situation, the exact configuration referring to some choice 
of the one-link bond weights considered below, some of which may vanish. 

A si bond particle of real-space coordinate fj involves A^^i two-site one-link bonds 
consistently with each family having Nsi/2D links of different type. The link type is 
labeled by an index (7 = 0, 1, ..., [Nsi/2D — 1] uniquely defined below. Each link of a si 
bond particle of real-space coordinate fj involves two sites of coordinates f — f^i and 
f+f^i where fj = f— f^i so that the link centre f^= fj + f^i is the middle point located 
half-way between the two sites and the link vector rj^ is defined below. The real-space 
coordinates fj = f— f^i and r + r^j belong to the sub-lattice 1 and sub-lattice 2 of the 
spin effective lattice, respectively. For each family there are Nsi/2D link vectors r|j 
which for the square lattice read, 

rli = rl + fl ; f°, = / 1 e., ; ^ = 0, 1, ..., [N,^/2D - 1] , (13) 

where d = 1,2, I = ±1, and T/^ is a T vector. It has Cartesian components 
T^i = [T^ii,T^i2] fo^^ the square lattice and Tf^ = [T-^^n] for ID. There are Nsi/2D 
T vectors T/^, one for each choice of the following Cartesian components, 

Tli,i = ^(^sNi- 1 = 1,2, 

Na =^,l,...,NajA-l- N^=-NajA + l,...,-l,^,l,...,NajA. (14) 

Here d = 1, 2, 1 = 2, 2 = 1, / = ±1 and A'",^ and A'j are consecutive integer numbers. 
The expressions provided in Eq. (fT3l) apply to the ID lattice as well provided that only 
the index value d = 1 is considered. The ID component of the T vectors is given by 
Tfi , = lasm where N, = 0, 1, ..., KJA - 1. 
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The link-type index ^ = 0, 1, ..., [Nsi/2D - 1] labels the Nsi/2D T vectors fj^^. For 
the square lattice it is defined in terms of the numbers N^ and N^ given in Eq. flT4|) and 
reads, 

g = Nd + 2\N^\^- iV,-<0, 

= iV, + 2(iV,--l)^; iV,->0. (15) 

For ID one has that g = Ni = Q,l, ..., A^,i/2 - 1. 

The values of the link- type index g are consecutive positive integers whose minimum 
value g = Q corresponds to A^i = A^2 = so that, 

Til = . (16) 

For the model on the square lattice the maximum value is (? = [Nsi/2D — 1] and refers 
to Nd = NajA - 1 and iVj = N^JA. 

Each pair of values (and value) of the Cartesian coordinates of the T vector 
^di = [^/zi!^/z2] for the square lattice (and T/^ = [T/^ ^] for ID) corresponds to 
exactly one of the values (? = 0, 1, ..., [Nsi/2D — 1] so that, 

E E - E i D = ^-: E - E ; « = i-(i7) 

Nd=0 N^=-Nas/'i+i 9=0 Ni=0 3=0 

Two-site links with the same g value and d ^ d' and/or / ^ I' are equivalent links. 
Those are of the same type but belong to different families. Furthermore, T vectors T/^ 
and T^, j, with the same value of g and d ^ d' and/or / ^ I' are related as follows. 



Td,i — ^^' [^d,d' + 5j,d'0"x] T^'^i, ; a^ 



1 

1 



where Cx is the usual Pauli matrix. 

As described below in terms of suitable operators, a si bond particle of real-space 
coordinate rj is a superposition of A^^i two-site one-link bonds each being associated 
with a link vector r|p For each site of the spin effective lattice there is exactly one 
other site of the same lattice such that the link connecting the two sites has centre 
at fj + rj^i- Therefore, any link of the same family involves two sites of well-defined 
real- space coordinate rj + fj^i — f^i and "f^j +r^i+r^i, which do not contribute together 
to any other link of the same family. 

An important quantity is the distance between the two sites of a link which we 
call link length or two-site bond length. It is independent of the real-space coordinate 
f = fj -\- f^i of the link centre and is fully determined by the link vector f^i and thus 
depends on the link type associated with the index g only. For D = 2 and D = 1 it 
reads, 

e, = |2f|J = aV(l + 2Ar,)2 + (2Ar,-)2; ^^ = |2xf^J = a, (l + 2Ari) , (19) 

respectively. Its minimum and maximum values are, 

mm^g = ^o = as-, max^g = V2a,{Naj2 ~ 1) + 0{1/Na) , (20) 



Ground-state two-spinon bonds in the Hubbard model on a square lattice 11 

for the square lattice and min^g = ^o = ds and niax^^ = a^ {Naj2 — 1) for ID. 

For ID, links with different g have different link length ^g. In turn, for the square 
lattice there are links of different type and hence different g which have the same link 
length ^g. Indeed, analysis of the link-length expression of Eq. ( !T9|l reveals that links 
with different g value and numbers [Nd, iVj] and [A^^, iVj], respectively, such that Nd = N'^ 
and Nd = —N'^ have the same link length. 

The set of values of the numbers A^^i and N2 given in Eq. flT^ imply that the 
maximum value of the link length is y/D as {Naj2 — 1) rather than \fD as {Na^ — 1). 
Indeed and as mentioned above, the links contributing to a si bond particle of the 
A^^^ = configuration state are independent. It is then required that each link involves 
two sites that participate simultaneously in exactly one of such links. Within the torus 
row and column periodic boundary conditions for the square spin effective lattice such a 
requirement is fulfilled provided that the range of the numbers A^^i and N2 is that given 
in Eq. (JH]). 

For each family of two-site links associated with a si bond particle of real-space 
coordinate r^ there is a primary link. It corresponds to g = and thus connects two 
nearest-neighboring sites of the spin effective lattice, one of them having the same real- 
space coordinate fj = r — r°; as the si bond particle. For primary links the link vector 



r 



I reads r^i = f^i where the primary link vector f^i is given in Eq. (TT3|) . 



For the model on the square lattice there are four primary links, one per family. 
Their link vectors r°; have components such that A^^i = A^2 = in Eqs. flT3|) and flT^ . 
Therefore, the primary links have minimum length f^o = a^. Alike the remaining links 

d.l 

of its family, the centre of a primary link is located at r = fj + f^i- For the square 
lattice there are two horizontal primary links whose centers are located at fj + rf; with 
/ = ±1 and two vertical primary links whose centers are located at fj + f2i with / = ±1. 
In the case of the ID lattice there are two primary links whose centers are located at 
fj + ffi with / = ±1. 

3.2. Partitions and g -primary partitions 

The building blocks of the A^j\ = configuration state are singlet pairs of spinous on 
sites f~ and f^ of the spin effective lattice, 

rj = r, + fl, T rl ; d = d{j) , / = l{j) , g = g{j) , (21) 

where the values of the integer indexes d, I, and g are in the ranges d = 1,2, I = ±1, 
and g G (0, Nsi/2D — 1), respectively, and are a function of the index j = 1, ..., Nsi of 
the real-space coordinate fj in the sublattice 1 of each si bond particle. Indeed, the two 
sites of such pairs of sites are connected by one-link bonds and each bond is associated 
with exactly one si bond particle. 

Each connection involving A^^i different bonds determines a partition. A partition 
is a 2A^si-spinon occupancy configuration where each site of the spin effective lattice is 
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Figure 1. Sub-domain of tlie spin effective lattice with a primary partition oi d = 1 
horizontal links for the model on the square lattice. The primary links are represented 
by the thick horizontal lines connecting two sites of the spin effective lattice. For a 
reference frame where the site located at the corner on the left-hand side and lower 
limit of the squared sub-domain has Cartesian coordinates (0,0) the family indices 
read d = I, I = +1 ii that site belongs to sub-lattice 2, whereas d ~ 1, I = —1 ii 
instead it belongs to sub-lattice 1. 



linked to one site only and all 2Nsi sites then correspond to A^^i well-defined two-site 
one-link bonds, each belonging to a different si bond particle. 

We recall that a si bond particle involves the superposition of Ngi such two-site 
one-link bonds whereas a partition involves one two-site one-link bond from each of the 
A^si si bond particles. Each si bond particle contributes with exactly one of its two-site 
bonds to a partition. In a partition any site of the spin effective lattice participates in 
one bond only and there is a single link attached to each site which connects it to some 
other site. And the latter site is attached to the former site only. 

The Ngi = configuration state is then represented as, 

A^.i Nsi 

Cpl[\r7,f;}; Cp = Ylh;,,, (22) 



p 



i=i 



i=i 



where the coefficients hg are associated with the bond weights and appear in the bond 



operators defined below, the product of singlet states Yl 



N,, 



T . 



') refers to a bond 



state associated with a given partition, and the summation Yip is over all partitions. 

A particular type of partition involves A^^i identical links. The indexes d, I, and g 
of identical links have the same values but correspond to si bond particles with different 
real-space coordinates fj. Such a partition involves a set of Ngi identical two-site one- 
link bonds whose links connect different sites of the spin effective lattice, each site being 
linked to exactly one site. In this case the two real-space coordinates of the A^^i pairs 
of sites are connected by the same real-space vector 2r|j so that each bond link has the 
same length. 

Each of the A^^i bonds of a partition involves two sites of real-space coordinates fj/ 
and fj/ + 'ifdi which belong to different sub-lattices where j' = 1, ..., A^si. The relation 
to the notation used above for the real-space coordinates of the two sites of a link is as 
follows. 



rr 



Tj + r 



dl 



dl I 



fji + 2r 



dl 



dl 



_l_ ^ 



dl ' 



J,J 



l,...,Ar,i, (23) 
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Figure 2. Sub-domain of tlie spin effective lattice with a primary partition of d = 2 
vertical links for the model on the square lattice. The primary links are represented by 
the thick vertical lines connecting two sites of the spin effective lattice. As in Fig. 1, 
the family indices read d = 2, / = +1 if the site of Cartesian coordinates (0, 0) belongs 
to sub-lattice 2, whereas d = 2, / = — 1 if instead it belongs to sub-lattice 1. 



where both fji and fj are real-space coordinates of the sub-lattice 1 and thus of the si 
effective lattice. Except for a primary link one has that j ^ j'. The site of real-space 
coordinate fj is that of the corresponding si bond particle. It is closest to the link 
centre at fj + f^i. In turn, fjt is the real-space coordinate of one of the two sites of the 
spin effective lattice involved in the link. 

When a partition is a set of A^^i identical primary links, all site pairs involve nearest- 
neighboring sites of the spin effective lattice. It is then called a primary partition. The 
family of a primary partition is labeled by the indexes d and / of the corresponding 
identical links. Figures 1 and 2 represent primary partitions oi d = 1 horizontal and 
d = 2 vertical links, respectively, for a sub-domain of the spin effective lattice of the 
model on the square lattice. 

An useful concept is that of a g -primary partition. It is defined as the superposition 
of the 2D = 2,4 primary partitions. It follows that a (yf-primary partition contains 
2D Nsi two-site primary links. In such a configuration each of the 2Nsi sites of the 
spin-effective lattice has 2D = 2,4 links attached to it. Figure 3 shows a sub-domain of 
the spin-effective lattice with the (^-primary partition of the N^^ = configuration state 
for the model on the square lattice. 



3.3. The si bond-particle operators and the subspace they act onto 

The spinon occupancy configurations considered above are similar to those associated 
with multi-spin wave functions of spin-singlet states used by several authors |lHl |20l |2l] ■ 
Such wave functions are often constructed having as building blocks two-site and two- 
spin spin-singlet configurations similar to that of Eq. (I2T]) except that here the two spins 
refer to sites singly occupied by rotated electrons and thus correspond to U /At > rather 
than only to U/4t ^ 1. In such schemes one also connects pairs of lattice sites with 
bonds and each such a connection determines a partition. However, here bonds involve 
sites of the spin effective lattice whereas those of previous related studies refer to the 
sites of the original lattice of the corresponding quantum problems. 

For a given partition one can define a valence bond state [HI [20|, [21] as a product of 
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Figure 3. Sub-domain of the spin-effective lattice representing tlie (/-primary partition 
of tlie N^i = configuration state for the model on the square lattice. The sites 
belonging to the sub-lattice 1 are represented by filled circles. The horizontal (and 
vertical) thin and thick lines refer to d = 1 (and d = 2) I — +1 and I = —1 primary 
links, respectively. For the square lattice each site of the spin effective lattice has in a 
(/-primary partition four links attached to it. 



singlet states and represent an arbitrary singlet by a superposition of valence bond states 
of general form similar to that of Eq. ( l22l) . That involves a sum over all partitions of the 
lattice into set of pairs. However, for general wave functions such a decomposition works 
in general very badly. Indeed, valence-bond states are not orthogonal and their basis is 
overcomplete. Fortunately, here each of the A^^i two-site bonds of a partition belongs 
to a different si bond particle so that each of such particles contributes to a partition 
with exactly one bond. Such a restriction eliminates the unwanted and unphysical 
contributions and renders the bond states free of the overcompleteness problem. 

For the one- and two-electron subspace and Na ^ oo, the operators grj,si (and 
9r si) which annihilate (and create) a si bond particle at a site of the spin effective 
lattice of real-space coordinate fj have the following general form both for the ID and 
square lattices, 

Nsi/2D-1 
9=0 



9f,,si 



D 



fj,sl,g — 2^ Z^ bf^j+f0^,sl,d,l,9 ' -O — I5 2 , 



(24) 



d=i i=±i 
so that the expression of gt ^j^ involves the operators at 



vt 



rj,sl,9 

y.gig and arj,si,g create and annihilate, respectively, a 



^t 



and h 



'f,sl,d,l,g 



are two-site 



(a,-^,,i,g)^ and6t,i,^;_ 
{brMAi^gV- The operators at 

superposition of 2D = 2,4 bonds of the same type and b^^^^, 
one-bond operators whose expression is given below. 

For the square (and ID) lattice the four (and two) primary links associated with the 
operators at ^-^ ^ and a^.^si.o are behind most of the spectral weight of a si bond particle 
of real-space coordinate Vj. Consistently, the absolute value \hg\ of the coefficients hg 
appearing in the expressions of such operators given in Eq. (^^ decreases for increasing 
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link length C,g- These coefficients obey the normalization sum-rule, 



[Nsi/2D-l] 

|2 -'- 



E 



^^' 2D' ^ = ^'2- (25) 

9=0 

The exact dependence of |/ig| on the link length ^g, value of U/At, and hole concentration 
X remains for the Hubbard model an involved open problem. The suitable use of this 
sum-rule and related symmetries leads though to useful information, as discussed below. 
By definition the real-space coordinates r,- of the si bond-particle operators of Eq. 1^^ 
are those of the sub-lattice 1. 

The two-site one-bond operators bl^-^^^, and bf=-si,d,i,g appearing in Eq. ( 1M|) are 
associated with a well-defined link connecting the two sites of real-space coordinates 
^ — ^d,i and r + r^^, respectively. Their expression can be obtained by considering the 
following related operator. 



-1 



st ;;t 



— -^[(1 - ^r-r|,,Jc._.s^^^C.^.,^^^ (1 - ^r+r|,,t 



— Jr~r^,,c Jr+r^,,c "r,sl,d,l,g ' (26) 

where in the second expression the operator b\,sidig reads. 



b^ 

r,sl,d,l,g 



:-iY 



V2 

t 



'1 

- + si -.s 
2 '■"'■'i,; 



r+r" 



d,l 



'\s^ 



(27) 



and bf^si,d,i,g = [brsidig] ■ Here the spinon operators are those given in Eq. ([6]). The 
second expression of Eq. ( 126|) is obtained by the use of Eq. ([8]). 

The phase factor (—1)'^^^ which appears in the operator of Eq. ( l27l) is associated 
with the d-wave symmetry of the si bond-particle two-spinon pairing of the model on 
the square lattice. The introduction of such a phase-factor refers to a self-consistent 
procedure which follows from the d-wave symmetry of the spinon energy dispersion 
found in Ref. [2] for the si fermions. Such objects emerge from the si bond particles 
studied here through a suitable extended Jordan- Wigner transformation and their 
energy- dispersion d-wave symmetry arises naturally from symmetries beyond the form 
of the operators introduced in Eqs. fl26|) and f l27|) . 

According to the above discussions, the real-space coordinates r — r^^ and r + r^i 
involved in the operators of Eqs. ( 126|) and fl27|l correspond to two sites that belong to 
different sub-lattices of the spin effective lattice, f = Tj + f^^ is the link centre, and 
the primary link vector r°; and link vector f^i are given in Eqs. flT3|) and f lT4|) . In the 
configuration generated by the operator of Eq. fl26|) the two sites are singly-occupied by 
the rotated electrons associated with the operators appearing in the ffist expression of 
that equation. 
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The N^^ = configuration state fl22p can be written in terms of si bond-particle 
operators given in Eq. fl24l) as follows, 

Nsi 

where |0s; A''^) is the spin SU{2) vacuum with A''^ = 2Nsi independent +1/2 spinous on 
the right-hand side of Eq. (jlj). It corresponds to a fully polarized spin- up configuration. 
The subspace where the operators of Eqs. (1Ml) - (!27j) act onto is defined by imposing the 
equality of the general configuration states given in Eqs. fl22l) and fl28|) . respectively. 
Such an equality implies several restrictions on the transitions generated by the two-site 
one-bond operators (!27j) which are summarized in four corresponding rules for exclusion 
of unphysical and unwanted spin configurations given below. Before introducing such 
rules let us discuss several properties of the si bond-particle operators which follow from 
the algebra given in Eqs. flA.7P - flA.9|) of Appendix A for the basic spinon operators of 



the two-site one-bond operators expressions provided in Eq. ( 127|) . 

The si bond-particle operators g\, ^^^ and g^p.^si of Eq. (!24l) involve a sum of A^^i two- 
site one-bond operators of general form given in Eq. (|27l) with Nsi/2D of such operators 
per family. The number of unoccupied sites A^^^ of Eq. ([9]) refers to a subspace with 
constant number A^^i of si bond particles. In turn, the creation and annihilation of one 
si bond particle by application of these operators, respectively, onto the ground state 
involves a superposition of A^^i elementary processes, which do not conserve the number 
of these objects. Each such an elementary process is generated by an operator h].^^^, 
and 6r,si,d,i,c/5 respectively, whose expression is given in Eq. fj27|) . 

Within the present LWS representation, application of the rotated-electron 
operators of Eq. ( l26l) onto two unoccupied sites of the original lattice generates 
two virtual processes. The first process involves creation of two c fermions and two 
independent +1/2 spinous. The second process refers to creation of a spin-singlet two- 
site and two-spinon configuration upon annihilation of two independent +1/2 spinous. 
Indeed, from analysis of the expression provided in Eq. ( l27l) . one finds that application 
of the operator h].^^^^ onto the sites of real-space coordinates f— rj'^ and r + f^i gives 
zero except when those sites are both occupied by an independent +1/2 spinon. This 
is consistent with a "unoccupied site" referring to two sites of the spin effective lattice 
with real-space coordinates r — r^^ and r + r^^, respectively, which are occupied by 
independent +1/2 spinous in the initial configuration state. 

As confirmed in Section 4 and consistently with the studies of Ref. [2], for 
excitations involving transitions between configuration states where a si bond particle 
moves around in the spin effective lattice by elementary processes which conserve the 
numbers A^^i and A^^^^, an independent +1/2 spinon plays the role of an unoccupied site 
of the si and spin effective lattices. In contrast, in elementary processes involving the 
creation of a si bond particle the two annihilated independent +1/2 spinous play the role 
of a "unoccupied site", which becomes occupied in the final occupancy configuration. 

According to the operator expression provided in Eq. (1271) . upon acting onto the 
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Figure 4. Small sub-domain of two g-primary partitions associated with Nj^^ = and 
Nj}^ = 2 configuration states for the model on the square lattice. In the configuration 
of the figure on the left-hand side the filled circle at the middle of the sub-domain 
corresponds to the site whose real-space coordinate fj is that of the si bond particle. 
Annihilation of that object leads to the configuration on the figure right-hand side. 
Note that for a primary g-basic partition this is equivalent to the suppression of the 
four links attached to the above site of real-space coordinate fj . 



independent-spinon occupancies the operator bl^^^, generates a superposition of two 
configurations. For one of these configurations the elementary process generated by that 
operator flips the spin of the spinon at site r + f^i and checks that the spin of the spinon 
at site f — f^i remains up. The elementary process generating the other configuration 
flips the spin of the spinon at site r — f^i and checks that the spin of the spinon at site 
f-\- r^i remains up. The relative phase factor —1 of the two configurations insures that 
the si bond particle created by the operator ^fl ^^ of Eq. ( !24l) is a suitable superposition 
of spin-singlet configurations. 

Figure 4 shows a small sub-domain of two (^-primary partitions of A^^^^ = and 
N^-^ = 2 configuration states. In the configuration on the figure left-hand side the 
filled circle at the middle of the sub-domain corresponds to a site whose real-space 
coordinate fj is that of a si bond particle. Application of the annihilation operator 
gi^^si of Eq. f l24|) onto that g-primarj partition occurs through the operator a^^^i^o also 
given in that equation. That leads to the configuration on the figure right-hand side. 
Hence annihilation of the si bond particle is for a gf-primary partition equivalent to the 
suppression of the four links attached to the above site of real-space coordinate fj. 

A superficial analysis of the configurations shown in Fig. 4 seems to indicate that 
there is one unoccupied site in the sub-domain of the final N!!-^ = 2 configuration state. 
However, if instead of the (/-primary partition one considers the corresponding four 
primary partitions one finds that there are two nearest-neighboring unoccupied sites. 
For the d = 1 and d = 2 primary partitions these two sites belong to the same row and 
column, respectively, as further discussed in Section 4. 

Concerning the application onto spin configurations of two-site one-bond operators 
and four-site two-bond operators, the restrictions arising from imposing that the 
representations of the N^-^ = configuration state given in Eqs. fl22|) and f l28|) . 
respectively, are identical correspond to the following four rules whose fulfillment 
prevents the generation of unwanted and unphysical spin configurations: 

First rule according to which application onto a spin configuration of a two-site 
one-bond operator bl^^^, or bfsia^i^g generates a physical spin configuration provided 
that its sites of real-space coordinates f—f^i and f+f^i, respectively, are in the initial 
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spin configuration (i) occupied by independent +1/2 spinous or (ii) linked by a bond. 

Second rule states that application onto a spin configuration of any of the elementary 
four-site two-bond operators &I-si,dJ,g^S'',si,d',«',(," hsi,d,i,gbl',si,d',i',g'^ hsi,d,i,gbf',si,d',i',g'^ 
and bt^^^i bf^si,d',i',g' gives zero when one of the two sites of real-space coordinates 
r — rj'j and r + f^i, respectively, is the same as one of the two sites of real-space 
coordinates r' — r^, j, and r' + f^, i,. Hence only when the two two-site one-bond operators 
of such four-site operators do not join sites or join both sites their application onto a 
spin configuration generates physical spin configurations. 

Third rule according to which when the two two-site one-bond operators 
of an elementary four-site two-bond operator 4,si,d,«,g4',si,d',r,g" b^,s'i-Ai,9bl',si,d',i',g'y 
br,si,d,Lgbf , si, d', I', g', and bl ^-^ ^i gbf ,si,d' ,1' ,g' ^ correspoud (i) to different si bond particles so 
that in the real-space coordinates f = fj + f^i and r' = fji + f^, ;, of their link centres, 
respectively, one has that j = j' or (ii) to the same si bond particle and to the same 
sites so that r ± r^^ = f ± f^, ;, the first rule applies independently to each of such 
two-site one-bond operators provided that in case (i) the second rule is obeyed. 

Forth rule refers to when the real-space coordinates of the link centres of the two-site 
one-bond operators of any of the four elementary four-site two-bond operators considered 
in the third rule are given by r = fj + f^i and f = fj + r °, ;,, respectively, so that 
such two-site one-bond operators correspond to the same si bond particle of real-space 
coordinate fj but f^^ ^ f^, ;, and states that then there are two cases. When in the 
initial spin configuration where the operators bi^^^^igbl,,.^^^,^,^, and bi^^si,d,i,gbl,^gi^d',i',g' 
(and bf^^si,d,i,gbf',si,d',i',g' and bi^^^^igbf^^si,d',i',g') act onto the sites of the spin effective 
lattice of real-space coordinates fj + r °, ;, — r J ;, and fj + r °, ;, + r J ;, are linked by 
a bond (and occupied by two independent +1/2 spinous) then a physical final spin 
configuration is generated provided that the sites of real-space coordinates fj +r^; ~'f^di 
and fj + r*°; + r'J'^ are (i) occupied by two independent +1/2 spinous or (ii) linked 
by a bond. In turn, when in the initial spin configuration where these four-site two- 
bond operators act onto the sites of the spin effective lattice of real-space coordinates 
fj + f^, II — f^i I, and fj + r^, ;, + f^, ;, are occupied by two independent +1/2 spinous 
(and linked by a bond) then a physical final spin configuration is generated provided 
that their sites of real-space coordinates fj + f^i ~ ^di ^"^^ ^i ~^ ^di ~^ ^di ^^^ linked by 
a bond (and occupied by two independent +1/2 spinous). 

Such rules follow naturally from the definition of the subspace where the operators of 
Eqs. ( l24l) -(!27l) act onto. The main criterion is that such operators have been constructed 
to inherently generating a faithful representation provided that the corresponding state 
( 128|) represents the N^^ = configuration state and hence is identical to that given 
in Eq. ( I22l) with the same value of A^^i = N^ . For instance, the second rule results 
from all partitions of the summation on the right-hand side of Eq. f l22|) a site of the 
spin effective lattice being linked to exactly only one site. Indeed, in a given partition 
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no two-site one bonds join the same site. Furthermore, the third rule refers to four- 
site two-bond operators whose two two-site one-bond operators belong to the same si 
bond-particle operator gjf.gi or (7^,51 yet correspond to different pairs of sites of the 
spin effective lattice. Such four-site two-bond operators appear in the expressions of 
the si bond-particle operators [gt ^-^j'^ or [5't=-_,-,si]^, respectively, which as confirmed in 
Appendix A give zero when acting onto the subspace which the operators of Eqs. fl2^ - 
fl27|) refer to. The point is that when in the initial spin configuration where the operators 

^,314,1,3^7^,314', I', g' ^^^ ^r,3l,d,l,gb^,^^^,i,g, (and bf^sl,d,l,gbf,3l,d',l',g' and b^^,i^d,l,g^r',3l,d',l',g') 

act onto the sites of the spin effective lattice of real-space coordinates Vj + r °, ;, 



■9 



V ;' 



] ^ ' d',l' ' d' 

and rj + r °, ;, + r|, j, are occupied by two independent +1/2 spinous (and linked by a 
bond) the operator on the right-hand side of S'^.^siS'^ ,si ^^^ gfj,sigl.^si i^^^ gfj,sigfj,si 
and gt gigfj,3i) which such four-site two-bond operators belong to, respectively, has 
transformed that configuration into an intermediate virtual state where such sites are 
linked by a bond (and occupied by two independent +1/2 spinous). 

The A^^^ = configuration state is given exactly by the same superposition of 
partitions when the real-space coordinates fj of the si bond-particle operators of Eq. 
fl24|) where j = 1,...,A'^^ are chosen to refer to any of the two sub-lattices of the 
spin effective lattice. That property is straightforwardly confirmed for the gf-primary 
partition represented in Fig. 3 and is fulfilled provided that one sums over all possible 
partitions. Hence, whether the above real-space coordinates are those of one or the 
other of such sub-lattices, the occupancy configuration associated with the creation 
of the Ngi = N^J2 bond particles is exactly the same. This confirms that the two 
corresponding choices for the si effective lattice describe the same A^^'\ = configuration 
state so that there is a gauge structure. This holds in the limit A^ ^ 1 that the present 
description refers to. 

3.4. Two-site one-bond weights 

The exact dependence of the absolute value \hg\ of the coefficients appearing in Eqs. 
fl22|) and fl25l) on the link length ^g remains an open problem. However, if one assumes 
that \hg\ has the following simple power-law dependence on that length, 

\hg\ = -^.. (29) 

the link-length expression (iT9l) and normalization condition (1251) alone imply that C^ 
be given by, 

^' ^ 4 Y.N, Eiv,[(l + 2Ar,)2 + (2Ar,)2]-.i ' ^ = 2' 

Here the summations are over the range of Nd and Aj (and A^i) values given in Eq. 
f fT7|) for the square (and ID) lattice. An expression of the general form (|29|) would 
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Figure 5. The relative spectral weight of the primary links R2D of Eq. (l32t for the 
square lattice as a function of the exponent asi for the range Ofgi > 1. That relative 
weight refers to the simple power-law expressions given in Eqs. 1^^ and ([50]) for the 
coefficients \hg\ appearing in the summations of Eq. (|24|) . The physical range of the 
exponent agi corresponds typically to i?2D > 0.9 so that the primary links are behind 
most of the si bond-particle spectral weight. For asi ~ 2.5 the ratio R2D is larger 
than 0.9 and approaches quickly the unit upon further increasing a^i. 



imply that the dependence of \hg\ on U/At and x occurs through that of the exponent 
tt^i = asi{U/4:t,x) only. 

For the coefficients \hg\ power-law expressions given in Eqs. ( I29l) and ( l30l) the ratio 
of the spectral weight of a primary link over the weight of the corresponding link family, 
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The ratios of Eq. (l32l) refer to the limit A^q — )■ 00 and are plotted in Figs. 5 and 6 for 
the square and ID lattice, respectively, as a function of the exponent a^i for the range 
ttsi > 1- For ttsi slightly larger than 5D/A where D = 1,2 the ratios Rin {D = 1) and 
R2D {D = 2) are larger than 0.9. 

For Na -^ 00 the exact coefficients \hg\ are decreasing functions of the link length 
whose expressions most likely are not, at least for the whole link- length range, of the 
simple form given in Eqs. (I29l) and (1301) . However the exact sum-rules (l25l) together 
with the coefficients \hg\ being decreasing functions of the link length implies that at 
least for link- length not too small the exact coefficients \hg\ fall off as in Eq. ( l29l) and 
the corresponding ratios of Eq. (!3T!) have a behaviour similar to that shown in Figs. 5 
and 6. Hence for the square lattice one expects that a^i > 2.5 so that the ratio R2D is 
larger than 0.9 and the primary links are behind most of the spectral weight of the si 
bond particle, alike for the ID lattice for a^i > 1.25. 

Here we considered the A^^'\ = configuration state where all sites of the si effective 
lattice are occupied. In the following we generalize our results to configuration states 
with a finite number of unoccupied sites. 
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Figure 6. The relative spectral weight of the primary links Rid of Eq. ((32)) for 
the ID lattice as a function of the exponent a^i for the range asi > 1. The ratio 
Rid approaches the unit upon increasing the value of a^i more quickly than the ratio 
plotted in Fig. 5 for the square lattice. For asi — Ugi — 1, 2, 3, ... integer the relative 
weight reads Rm = 1/[(1 — 2^^"^0C(27^si)] where C{x) is a Riemann's zeta function. 
For instance Rid = S/vr^ w 0.811 for q^i = 1 and Rid = 96/tt'^ « 0.986 for a^i = 2. 



4. The N^-^ = 1,2 configuration states and the kink-hke and anti-kink-hke 
link occupancies associated with the unoccupied sites 

The one- and two-electron subspace considered in this paper is spanned by states with 
none, one, and two si effective-lattice unoccupied sites. There are several configuration 
states with one or two unoccupied sites in that lattice so that such sites can move around 
in it. We recall that such a motion is independent of that of the rotated electrons that 
singly occupy sites of the original lattice relative to the rotated-electron unoccupied and 
doubly-occupied sites of that lattice. The latter are instead described by the motion of 
the c fermions whose occupancy configurations correspond to the state representations 
of the charge global U{1) symmetry. In turn, the configuration states associated with 
the motion around in the si and effective lattices of the si bond particles relative to 
their unoccupied sites refer to state representations of the global spin SU{2) symmetry. 
The global r^-spin SU{2) symmetry does not play any role in the one- and two-electron 
subspace, since it corresponds to a single occupancy configuration (x > 0) or none 
occupancy configurations {x = 0) for all states which span that subspace. For x > 
this is the r^-spin vacuum on the right-hand side of Eq. (jl]), which is invariant under the 
electron - rotated-electron unitary transformation. For x = one has that N^ = so 
that the ?7-spin effective lattice does not exist. 

Analysis of the transformation laws of the objects whose occupancy configurations 
generate the energy eigenstates of the square-lattice model in the one- and two-electron 
subspace p] reveals that the excited states generated by application onto a; > and 
m = ground states of one- and two-electron operators have one unoccupied site 
and none or two unoccupied sites in the si effective lattice, respectively. The energy 
eigenstates that span such a subspace are generated in Ref. [2] by c and si fermions 
momentum occupancy configurations. The si fermion occupancies in the si momentum 
band introduced in Ref. [2j that generate energy eigenstates with one and two si band 
holes can be expressed as suitable superpositions of the set of iVJ\ = 1 and N^-^ = 2 
configuration states, respectively, studied in the following. 
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The structure of the N^^ = configuration states is simpler than that of the 
N!!^ = 1,2 configuration states. A key property simphfies the study of the latter states: 
There is an one-to-one correspondence between the occupancies of the Ngi occupied sites 
of the si effective lattice of the partitions of a iVJ\ = configuration state and those of 
the occupied sites of the si effective lattice of partitions of the N^-^ = 1,2 configuration 
states with the same number Ngi of si bond particles. The point is that except for a 
change in the real-space coordinates of some of the sites, the link configurations of the 
sites of the spin effective lattice in one-to-one correspondence with each other remain 
unaltered. 

4-1- A^s'\ = 1 g-primary partitions of the four sub-configuration states 

We start our analysis by considering the configuration states with one unoccupied site. 
The si fermion momentum occupancy of excited [A^-^ — N^] = 1 states of x > and 
m = ground states generated by application of one-electron operators onto the latter 
states is described by a suitable superposition of the set of N^-^ = 1 configuration states 
studied in the following. For such configuration states the unoccupied site is associated 
with their single independent +1/2 spinon. For simplicity we have considered above 
that for the N^^ = configuration state the number of sites A^^ = 2Sc = Nc of the 
spin effective lattice is such that A'^^ is ior D = 2 an integer so that the spin effective 
lattice is a square lattice. In turn, for configuration states with one unoccupied site the 
number A^^ of sites of the spin effective lattice is odd so that for D = 2 the designation 
A"^ is not to be understood to imply that the number A^^^ is an integer. 

We recall that the number of sites of the spin effective lattice is given by A"^ = 
2Sc = Nc where Sc is the eigenvalue of the generator of the global U{1) symmetry 
Sc = V^ ScV, Sc = Q/2, the operator Q is given in Eq. ([3]), and Nc is the number of 
c fermions. Since Sc is a half-off-integer, according to the studies of Ref. [6] the global 
S0{3) X 5*0(3) X t/(l) symmetry of the present quantum problem implies that the r^-spin 
Srj and spin Ss are half-off-integers as well. It follows that the number of sites of the si 
effective lattice is according to Eq. (|9]) given by A"^^ = [A'i^/2 + Ss] where Sg = 1/2. 
Out of those, A^si = [A'i^j — Ss] are occupied by si bond particles and one is unoccupied, 
consistently with A^^^^ = 2Ss = 1. Hence such states have also one unoccupied site in 
the si effective lattice. 

For the study of the A^j\ = configuration state we considered in Section 3 a 
change of gauge structure such that the real-space coordinates of the sites of the si 
effective lattice correspond to one of the two sub-lattices of the spin effective lattice. 
The results of Ref. [2] confirm that for the limit N^ ^ 1 and hole concentrations x such 
that (1 — x) > l/N^ the two choices of si effective lattice lead to the same quantum 
numbers for the whole one- and two-electron subspace. 

Let us consider A^^^^ = 1 configuration states whose spin effective lattice has 2 A^^i + 1 
sites and for the model on the square lattice the spin effective lattice of the corresponding 
N^i = configuration state with the same number A^^i of si bond particles is a square 
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lattice and has Na^ x Na^ = 2Nsi sites. For the present A^^ — > oo hmit the correct 
physics is then achieved if one considers that out of the 2Nsi + 1 sites of the spin 
effective lattice of such A^^^^ = 1 configuration states, 2Nsi sites correspond to a square 
lattice and the position of the extra site is well defined and given below. The latter site 
does not belong to the square lattice formed by the remaining 2Nsi sites of the spin 
effective lattice and has suitable boundary conditions such that it is the unoccupied site 
of the si effective lattice. However, for the spin effective lattice the extra site and the 
unoccupied site are not always the same site, as discussed below. 

For the square (and ID) lattice the extra site boundary conditions are compatible 
and consistent with the periodic boundary conditions of each row and column (and 
chain). For the A^^^ = configuration state of the model on the square lattice such a 
row and column periodic boundary conditions imply torus periodic boundary conditions 
for the spin effective lattice. Here the presence of the extra site slightly affects the latter 
boundary conditions. Hence the spin effective lattice of iVJ\ = 1 configuration states has 
no pure torus periodic boundary conditions. Within the extra-site boundary conditions 
given below each N!!^ = 1 configuration state has for the square (and ID) lattice four 
(and two) sub- configuration states. The spin effective lattice is the same for all such 
sub-configuration states. The concepts of a partition and (yf-partition remain the same 
as for the corresponding N^-^ = configuration state, except that now there is an extra 
site. 

For the model on the square lattice the extra site of the spin effective lattice belongs 
to the row or column that for a. d = 1 and d = 2 primary partition specified below, 
respectively, the unoccupied site belongs to. The real-space coordinate of the extra site 
and its Cartesian coordinates are denoted by. 

In the present Na —^ oo limit the following periodic conditions yield the correct physics 
and assure that the extra site is uniquely defined and the same site for all partitions of 
a given A^^^ = 1 configuration state, 

a;° = x° + /(L + a); 4 = x°-; rf=l,2,/ = ±l, (34) 

where the index rf = 1, 2 is defined as in Eq. JHM . Such boundary conditions apply to 
the square lattice for d = 1,2 and / = ±1 and to the ID lattice for d = 1 and Z = ±1. 
For the row {d = 1) [and column {d = 2) for the model on the square lattice] the extra 
site belongs to, they are equivalent to periodic boundary conditions. Indeed that row 
[and column] has A"a -|- 1 sites and hence its length is L -|- a = {Na + 1) a. 

For each sub-configuration state there is exactly one (/-primary partition, which is a 
superposition of the corresponding four primary partitions. Hence instead of 2D = 2,4 
primary partitions, for a A^j\ = 1 configuration state there are 8D = 8, 16 primary 
partitions. For all partitions of the summation on the right-hand side of Eq. (!22ll 
the extra site has the same Cartesian coordinates. Nonetheless, the extra site and the 
unoccupied site are not the same site of the spin effective lattice for all such partitions. 
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For instance, the extra site and the unoccupied site are the same site for three out of 
the four primary partitions of a primary gf-basic partition, as confirmed below. 

For the model on the square lattice the periodic boundary conditions of the sites 
belonging to rows and columns other than that of the extra site refer to a length L rather 
than L + a. For the extra site the boundary conditions ( 1M1) apply to both the square 
and ID lattice and are valid for the 2D = 2,4 sub- configuration states considered below. 
In turn, for the model on the square lattice and concerning the remaining sites of the 
extra-site row and column such periodic boundary conditions apply to those belonging 
to the row for horizontal sub-configuration states and to the column for vertical sub- 
configuration states, respectively. For horizontal (and vertical) sub- configuration states 
the sites of the extra-site column (and row) obey periodic boundary conditions associated 
instead with the length L, alike the sites that do not belong to the extra-site column and 
row. That the periodic boundary conditions of one row or column of the spin effective 
lattice refer to the length L + a rather than L is for the model on the square lattice 
behind the above deviation from the pure torus periodic boundary conditions. 

Except for the unoccupied site, the remaining 2A'^si sites of the spin effective lattice 
are occupied by A^^i si bond particles. Furthermore, out of the 2Nsi + 1 sites of 
that lattice, for the 2Nsi sites other than the extra site, which for some partitions 
is not the unoccupied site, there are two well-defined sub-lattices, alike for the A^^^ = 
configuration state. Indeed, the extra site belongs to the spin effective lattice but does 
not belong to any of these two sub-lattices. Moreover, now sub-lattice 1 is that whose 
real-space coordinates are the same as those of the sites of the si effective lattice whose 
row or column are not that of the extra site. In turn, for horizontal (and vertical) sub- 
configuration states the occupied sites of the si effective lattice of the extra-site row 
(and column) belong instead to sub-lattice 2. 

The si bond-particle operators act onto a subspace with constant values A^^ = 
2Sc = Nc of spin-effective-lattice sites. Therefore, suitable analysis of the general 
expressions of the number of sites and unoccupied sites of the au effective lattices 
[2] reveal that for au = si one of these operators or products of such operators can 
generate transitions where the value of N^^ either is conserved or changes by an even 
number. Indeed, transitions where such a value changes by an odd number imply that 
that of the number A*"^ = 2Sc = Nc of both spin-effective-lattice sites and c fermions 
also changes. For instance, transitions from the A^^^^ = configuration state to a A^^^ = 1 
configuration state involve creation or annihilation of one electron. Hence they involve 
creation or annihilation of one rotated electron as well. In terms of the objects of our 
description that process involves creation or annihilation, respectively, of a c fermion 
plus creation of an independent spinon [2]. In the LWS representation the latter is an 
independent -1-1/2 spinon whose creation is equivalent to adding an extra site to the 
spin effective lattice. 

Therefore, addition of the extra site is a process beyond the si bond-particle algebra 
whose operators are well defined for subspaces with constant values of A^^ = 2Sc = Nc 
spin-effective-lattice sites only. Fortunately, our general description also accounts for 
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such a process, which involves both creation or annihilation of a c fermion and creation 
of an independent +1/2 spinon. In turn, the generators of the transitions between 
different A^^'\ = 1 configuration states involve si bond-particle operators only. Such 
transitions are equivalent to the motion around in the spin effective lattice of si bond 
particles. But since in the present case there is a single unoccupied site it is often more 
convenient to describe such a motion in terms of those of that site. 

The unoccupied site has the same position in the si effective lattice and gf-primary 
partitions of the 2D = 2,4 sub- configuration states. Such a position is the same as 
that of the extra site in the spin effective lattice. For each position of the extra 
site in that lattice there is a A^^^^ = 1 configuration state. Specifically, there are 
^asi — Wasf^ + ^s] = ^si + 1 such Configuration states. The si fermion occupancies of 
Ref. [2] of the above mentioned excited states generated by application onto x > and 
m = ground states of one-electron operators can be expresed as suitable superpositions 
of such A^^^ = [N^J2 + Ss] = Nsi + 1 configuration states. For the model on the square 
lattice the four sub-configuration states of each of such A^^^ = 1 configuration states 
are called horizontal-kink, horizontal- anti-kink, vertical-kink, and vertical- anti-kink sub- 
configuration state. 

It useful for the description of the A^^^^ = 2 configuration states of Subsection 4-5 to 
introduce the A^j\ = 1 g-pnmaxj partition of the horizontal-kink sub- configuration state 
as a limiting case of a general configuration whose unoccupied site and extra site do not 
have the same real-space coordinate in the si effective lattice. The (/-primary partition of 
the above sub-configuration state is then reached by considering the particular case when 
the unoccupied and extra site are the same site. As discussed below, the unoccupied 
site of a ^f-primary partition (and that of the underlying si effective lattice) is identified 
with the unoccupied site of three out of its four primary partitions. A sub-domain of the 
spin-effective lattice with a part of such a general gf-primary partition is shown in Fig. 7 
for the square lattice. Its extra site is not shown in the sub-domain under consideration 
and is not the unoccupied site marked by the open circle. The occupied sites of the si 
effective lattice are represented by filled circles. 

The row kink-like link occupancy starting on the right-hand side of the unoccupied 
site ends at the extra site. We say that such a link occupancy and those of the 
surrounding rows have opposite signs. Indeed, the si bond-particle sites of the row part 
associated with the former occupancy belong to sublattice 2 of the spin effective lattice 
whereas those of the si bond particles on the surrounding rows belong to sublattice 1. 

The (7-primary partition of the horizontal-kink sub-configuration state corresponds 
to a particular case of that represented in Fig. 7 whose extra and unoccupied site are the 
same site so that the kink-like link occupancy extends over the whole unoccupied-site 
row. Figure 8 shows a sub-domain of the spin-effective lattice representing the ^f-primary 
partition of the N^^ = 1 horizontal-kink sub-configuration state for the square lattice. 
Since one of the issues to be further discussed is the periodic boundary conditions of the 
extra-site row, for simplicity the A^^^ ^ 1 sites of a row are represented in Fig. 8 by a 
few sites. The last site of each row is marked by a X to indicate that it is the same site 



Ground-state two-spinon bonds in the Hubbard model on a square lattice 



26 




Figure 7. Sub-domain of tlie spin-effective lattice representing a general type of link 
occupancy for the model on the square lattice which when the extra and unoccupied 
sites have the same real-space coordinate corresponds to the g-primary partition of the 
iVJ\ = 1 horizontal-kink sub-configuration state. However, in the link configuration 
shown here the unoccupied site does not coincide with the extra site so that the 
former site plays the role of a mobile domain wall, alike for some g-primary partitions 
of Nj^i = 2 configuration states considered in Subsection 4-5. The horizontal (and 
vertical) thin and thick lines refer to d = 1 (and d = 2) I — +\ and I = —1 primary 
links, respectively. For the d = 1, / = +1 partition the unoccupied site is that attached 
to three links. In turn, for the d ~ 2, I — ±1 partitions and d — 1, I ^ —1 partition 
the unoccupied site is that which has a single link attached to it. The latter site is 
marked by an open circle in the figure and plays the role of unoccupied site of the si 
effective lattice. The occupied sites of that lattice are represented by filled circles. The 
extra site does not belong to the sub-domain shown here and marks the end of the 
horizontal kink-like link occupancy located on the right-hand side of the unoccupied 
site. 



as the first site of the same row. In turn, for the vertical direction the figure represents 
a sub-domain only so that no equivalent sites are included, alike for the sub-domain 
shown in Fig. 7. Note that for the extra-site row of Fig. 8 the extra site is marked by a 
X and thus is the same as the unoccupied site marked by an open circle. Therefore and 
in contrast to Fig. 7, in Fig. 8 the kink-like link configuration extends over the whole 
extra-site row. 

Both for the d = 2, I = ±1 primary partitions and d = 1, I = —1 primary partition 
of the (^-primary partition shown in Fig. 8 and that of Fig. 7 the unoccupied site 
is the site attached to a single link, marked by the open circle in both figures. The 
d = 1, I = +1 primary partition is the only one for which the unoccupied site is that 
attached to three links in the figures. The site marked by the open circle plays the role 
of the unoccupied site of both the (yf-primary partition and si effective lattice of the four 
sub-configuration states. 

For the g-priraarj partition of the horizontal-kink sub-configuration state we call 
kink row or 1, +1 line the link occupancy of the sites belonging to the same row as 
the unoccupied site. The notation d, I line where d = 1 and I = +1 stems from the 
d = 1, I = +1 primary partition being the only one out of the four such configurations 
of the corresponding gf-primary partition for which the unoccupied site is attached to 
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Figure 8. Sub-domain of tlie spin-effective lattice representing the ^-primary partition 
of the NJJ'i — 1 horizontal-kink sub-configuration state for the model on the square 
lattice. The horizontal (and vertical) thin and thick lines refer to d = 1 (and d — 2) 
I = -1-1 and I = —1 primary links, respectively, the unoccupied site is marked by an 
open circle, and the sites corresponding to real-space coordinates of si bond particles 
by filled circles. For the vertical direction only part of the configuration is shown. In 
turn, for the horizontal direction the Na^ ^ 1 sites of a row are represented by a few 
sites. Indeed, one of the goals of the figure is to illustrate the row periodic boundary 
conditions, a site marked by a X being the same as the first site on the left-hand site of 
the same row. For instance, the unoccupied site marked by an open circle is the same 
site as the last site on the right-hand site of the same row. Note that the extra-site 
row has one more site than the remaining rows and such a site is both the extra site 
and the unoccupied site so that the kink-like link configuration extends over the whole 
extra-site row, in contrast to the configuration shown in Fig. 7. On the other hand, 
the link configurations associated with the unoccupied site are the same as in Fig. 7. 
For the d — 1, ^ = +1 primary partition the unoccupied site is that attached to three 
links and for the d = 2, I — ±1 primary partitions and d = 1, Z = — 1 primary partition 
that which has a single link attached to it. 



three links. For the g-phraaiy partition of each of the four sub- configuration states of 
a N^i = 1 configuration state there is a d,l hne and in each case the corresponding 
d, I primary partition is that for which the unoccupied site is attached to three hnks. 
Furthermore, for the remaining three primary partitions the unoccupied site is always 
that attached to a single link: It is the extra site, as illustrated in Fig. 8. 

We call d, I g-primary partition that which involves a d, I line. Let fj^^ denote 
the real-space coordinate of the unoccupied site in the spin effective lattice for each of 
the four primary partitions of the d, I gf-primary partition. For each of such primary 
partitions it is given in terms of the real-space coordinate fj^ of the extra site provided 
in Eq. ( l33l) as follows, 

r*,Q ; d, ±/ and d, —I primary partitions , 

fjg + 2r^i ; d, I primary partition , (35) 

where the index d is defined as in Eq. (IT4|) . 

Figure 9 shows four small sub-domains of the spin-effective lattice representing part 
of general ^f-primary partitions of the type represented in Fig. 7. As discussed below in 
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Figure 9. Four small sub-domains of the spin-effective lattice showing part of g- 
primary partitions. Alike in Fig. 7, each sub-domain includes an unoccupied site 
which plays the role of a mobile domain wall. The first two configurations of the figure 
refer to kink-like and anti-kink-like link configurations, respectively, of horizontal sub- 
configuration states. The two last configurations correspond to kink-like and anti-kink- 
like link configurations, respectively, of vertical sub-configuration states. Such mobile 
domain walls occur in some N^i = 2 sub-configuration states with two unoccupied sites 
on the same row or column rather than in N^i — 1 sub-configuration states. However, 
the link structure of the unoccupied site is the same as for the latter states. 



Subsection 4-5, such sub-domains involve unoccupied-site mobile domains walls which 
occur in N^^ = 2 sub-configuration states when the two unoccupied sites belong to the 
same row or column. However, the link structure associated with the unoccupied site at 
the center of each of such sub-domains is that also occurring in the (yf-primary partitions 
of the four iVJ\ = 1 sub-configuration states, respectively. In turn, each of these sub- 
domains includes four primary links of a si bond particle with a single link attached 
to the unoccupied site. As confirmed for the horizontal-kink sub-configuration state 
on comparing it with the link occupancy configurations of Fig. 8, the only difference 
relative to the A^^^^ = 1 sub-configuration states refers to the structure and position of 
the remaining three primary links of that si bond particle, the link-structure of the 
unoccupied site being the same as that shown in Fig. 9. 

Specifically, the unoccupied-site link structure of the top configurations on the left- 
and right-hand side of Fig. 9 are the same as those of the g-priraaxj partitions of A^^^^ = 1 
horizontal-kink and horizontal-anti-kink sub-configuration states, respectively. For the 
d = 2, I = ±1 primary partitions and d = 1, I = —1 (and d = 1, I = +1) primary 
partition on the left-hand side (and right-hand site) of the figure top the unoccupied 
site is that attached to a single link, marked by an open circle. In turn, for the d = 1, 
I = +1 (and d = 1, I = —1) primary partition the unoccupied site is that attached 
to three links. Moreover, for the d = 1, I = ±1 primary partitions and the d = 2, 
I = —1 (and d = 2, I = +1) primary partition on the left-hand side (and right-hand 
side) of the figure lower-limit associated with the vertical-kink (and vertical-anti-kink) 
sub- configuration state the unoccupied site is that attached to a single link, marked by 
an open circle. On the other hand, for the corresponding d = 2 and I = +1 (and d = 2, 
/ = — 1) primary partition the unoccupied site is that attached to three links. Such four 
types of unoccupied-site link structures correspond to the gf-primary partitions of the 
four N^i = 1 sub-configuration states, respectively. 



Ground-state two-spinon bonds in the Hubbard model on a square lattice 29 

The (yf-priinary partitions of the horizontal-kink and horizontal-anti-kink (and 
vertical-kink and vertical-anti-kink) sub-configuration states include a kink row as that 
represented in Fig. 8 and a anti-kink row (and a kink column and anti-kink column), 
respectively. Such a row (and column) link occupancy configuration involves the set of 
sites belonging to the same row (and column) as the extra site. As mentioned above, 
an alternative designation for the kink row, anti-kink row, kink column, and anti-kink 
column is 1, -|-1 line, 1, —1 line, 2, -|-1 line, and 2, —1 line, respectively. The general d, I 
line occurs in a primary (yf-basic partition for whose d, I primary partition the unoccupied 
site is attached to three links. 

There is an apparent contradiction concerning the following issue. On the one hand, 
that the extra-site periodic boundary conditions ( 1M|) refer to both a row and column 
of Na^ + 1 sites is valid for the 2D = 2,4 sub-configuration states simultaneously. 
On the other hand, concerning the (yf-primary partitions of such states, Na^ + 1-site 
periodic boundary conditions apply to either the sites of the spin effective lattice other 
than the extra site belonging to the extra-site row or extra-site column, respectively. 
Specifically, for gf-primary partitions of horizontal (and vertical) sub-configuration states 
such periodic boundary conditions apply to the sites belonging to the 1, / (and 2, /) line 
only, whereas for the sites belonging to the extra-site column (and row) other than the 
extra site the periodic boundary conditions refer to A''^^ sites rather than to A''^^ + 1 
sites. 

The reason behind that apparent inconsistency is that the link occupancy 
configurations of each of the 2D = 2,4 d,l lines of the 2D = 2, 4 (yf-primary partitions, 
respectively, are different and thus refer to different sub- configuration states yet they 
all involve the same extra site. We recall that the extra site of Figs. 7-9 is linked to one 
site only. For fyf-primary partitions with a 1, / (and 2, /) line the extra site is linked to a 
site belonging to that line only so that only the real-space coordinates of the remaining 
sites belonging to the 1, / (and 2, /) line are shifted, whereas those of the sites belonging 
to the extra-site column (and row) remain unaltered. 

In addition to the d, I line of a (/-primary partition, it is useful to consider its two 
nearest neighboring lines. For each such an occupancy configuration there are three 
d,l,l' lines where /' = — 1,0, +1. Here d,l,0 = d,l line is a more general notation for 
the d, I line. For a kink or anti-kink row the 1, /, /' line is the row above it for /' = +1 
and below it for /' = —1. Moreover, for a kink or anti-kink column the 2, /, /' line is the 
column on its right-hand side for /' = +1 and on its left-hand side for /' = — 1. 

4-2. Motion of the unoccupied site around in the spin effective lattice 

The motion of the unoccupied site around in the spin effective lattice discussed here can 
be described as well in terms of the motion of the corresponding unoccupied site of the 
si effective lattice around in such a lattice. In the case of the si fermion description of 
the square-lattice model investigated in Ref. [2] the generators of the motion considered 
here are the two si translation generators in the presence of the fictitious magnetic field 
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Bgi associated with the Jordan- Wigner transformation that maps the si bond particles 
into si fermions. The components Qj^^ and Qj^^ *^f ^^e discrete momentum values Qj 
of the si fermions are the eigenvalues of such two si translation generators. Since the 
si fermion creation and annihilation operators are defined in and act onto subspaces 
spanned by mutually neutral states [2] such operators commute and the two components 
qj^^ and qj^^ can be simultaneously specified. 

It follows from the above analysis that for the si effective lattice and the 2D = 2,4 
^f-primary partitions of all sub-configuration states the unoccupied site has a well-defined 
position. For A^j\ = 1 there are A^^^-^ = A^^i + 1 configuration states and this is also 
the number of different positions that the unoccupied site can have in the si effective 
lattice. 

Transitions between configuration states where the unoccupied site has different 
positions in the si effective lattice result from a set of well-defined elementary steps 
which are behind its motion around in the spin and si effective lattices. The generators 
of such processes have simple and obvious expressions in terms of si bond-particle 
creation and annihilation operators. The above elementary steps correspond to changes 
in the occupancy configurations of rotated electrons and the unitary operator V 
preserves the occurrence of nearest hopping only for such objects. Therefore, in the 
case of the square lattice rotated electrons can hop vertically or horizontally only and 
given the relation between the original lattice and the spin effective lattice that leads 
to restrictions in the changes of the position of the link-site attachments. Indeed, the 
elementary steps behind the net motion of the unoccupied site involve horizontal or 
vertical shifts by ±0^ of some of such attachments where a^ is the lattice constant of 
the spin effective lattice given in Eq. ([T]). The net length passed by the unoccupied site 
in an elementary step is 2as- 

For ID the unoccupied site moves by ±20^ in each elementary step. As for 
the square-lattice case, such elementary steps involve smaller virtual elementary steps 
where link-site attachments move by zha^. For the square lattice the net motion of the 
unoccupied site is given by ±2as e^^ where i = 1 and i = 2 for horizontal and vertical 
elementary processes, respectively. In such elementary steps the net length 2as passed 
by the si bond particle corresponds to horizontal or vertical motion. In addition, for the 
square lattice there are horizontal- vertical or vertical-horizontal elementary steps where 
the unoccupied site moves by ia^ [e^;- -|- Cx.,] or ia^ [e^;- — Cx.,] where i = 1,2 and i' 7^ i. 
For the latter elementary steps the net length 2as corresponds to both a horizontal and 
a vertical ia^ shift, so that the initial and final positions distance is \/2as 

For the gf-primary partitions any elementary process involves the virtual 
annihilation of one of the four si bond particles surrounding the extra site. That 
virtual process involves the creation of two unoccupied sites and is followed by suitable 
horizontal and/or vertical zba^ shifts of well-defined link-site attachments. Those shifts 
are equivalent to the unoccupied site moving by ±2as txi-, icts [ex^+Ca;^,], or ios [e^;. — e^;^,] 
where i = 1, 2 and i' 7^ i. Finally, a si bond particle is created onto the site of the si 
effective lattice unoccupied in the initial configuration so that two unoccupied sites of 
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the virtual state are occupied. 

For some g-phmaij partitions of the model on the square lattice such elementary 
steps also involve collective processes where a whole kink row or anti-kink row (and kink 
column or anti-kink column) configuration moves from one row (and column) to the 
nearest-neighboring or second- nearest-neighboring row (and column). Such processes 
consist of a collective horizontal (and vertical) shift of —Og or +0^ of all the sites of 
a kink or anti-kink row (and kink or anti-kink column), respectively, plus a collective 
horizontal (and vertical) opposite shift of +as or —Os of all the sites of a nearest- 
neighboring or second-nearest-neighboring kink or anti-kink row (and kink or anti-kink 
column), respectively. 

Specifically, for (7-primary partitions of horizontal (and vertical) sub-configuration 
states, unoccupied-site horizontal (and vertical) elementary steps do not involve such 
kink or anti-kink row interchanges (and kink or anti-kink column interchanges). 
In contrast, for (^-primary partitions of vertical (and horizontal) sub-configuration 
states, unoccupied-site horizontal (and vertical) elementary steps involve kink or anti- 
kink column interchanges (and kink or anti-kink row interchanges) of second-nearest- 
neighboring columns (and rows). On the other hand, horizontal- vertical or vertical- 
horizontal elementary steps involve kink or anti-kink row interchanges of nearest- 
neighboring rows for g-priraarj partitions of horizontal sub-configuration states and 
kink or anti-kink column interchanges of nearest-neighboring columns for primary g- 
primary partitions of vertical sub- configuration states. 

Any motion of the unoccupied site around in the spin effective lattice can be 
generated by a set of the above elementary steps. The processes behind such steps 
are easiest to describe for the primary partitions. However, they are also well-defined 
for the remaining partitions. 

4.3. Correspondence between the N^ = 1 and N^ = configuration states 

The partitions that span the 2D = 2,4 sub-configuration states of a N^^ = 1 
configuration state are independent so that the number of partitions of the latter state is 
2D = 2,4 times larger than that of the N^-^ = configuration state. Indeed, there is an 
one-to-one correspondence between the partitions of that state and those of a A^^^^ = 1 
sub-configuration state with the same number of si bond particles. Concerning the 
primary partitions of N^^ = 1 sub-configuration states, except for the sites of the spin 
effective lattice belonging to the row and column which the unoccupied site belongs to, 
the real-space coordinates of the sites of that lattice have the same real-space coordinates 
as those of the corresponding primary partitions of the A^j\ = configuration state with 
the same number of si bond particles. 

A Ngi = 1 configuration state is defined by the position of its extra site whose 
real-space coordinate in the si effective lattice is the same as that of the corresponding 
unoccupied site. Except for the sites belonging to the kink row and anti-kink row (and 
kink column and anti-kink column), respectively, for the horizontal-kink and horizontal- 
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anti-kink (and vertical-kink and vertical-anti-kink) sub-configuration state the two sites 
of each of the A^^i links of the 2D = 2,4 primary partitions have the same real-space 
coordinates as the corresponding primary partitions of the N^^ = configuration state. 
The shifts of the real-space coordinates of the sites of the primary partitions of the 
Ngi = configuration state which become sites oi a. d,l line in the corresponding final 
sub-configuration state lead to the following real-space coordinates for the latter sites, 

r!^ -^ r- = r^ + 2f °j = r,° + / a, q . (36) 

Here r? denotes the corresponding real-space coordinate of the initial A^j\ = 
configuration state. 

The procedure to construct the non-primary partitions of a A^^^ = 1 sub- 
configuration state profits from all two-site links remaining exactly the same and only 
the real-space coordinates of the sites of the spin effective lattice of the A^^^^ = 1 sub- 
configuration state whose real-space coordinates are shifted for the primary partitions 
as given in Eq. (!36|l being shifted exactly by the same amount for the corresponding 
non-primary partitions. Such a procedure applies to any of the A^^i + 1 positions of 
the unoccupied site in the si effective lattice and hence to all corresponding A^^^^ = 1 
configuration states. It corresponds to the elementary processes generated by application 
of an one-electron operator onto a x > and ?ti = ground state. 

4.4- si bond-particle operators of N^ = 1 configuration states 

Each of the 2D = 2,4 sub-configuration states of a A^^^^ = 1 configuration state has a 
d, I line in the corresponding d, I (/-primary partition where for ID the index d is given 
only hj d = 1. We use a single index i = d,l to denote the indices d, I which are used 
to label each sub-configuration state according to the type of d, I line in its (/-primary 
partition. For a A^^^^ = 1 configuration state the operators given in Eq. (^^ for the 
N^i = configuration state read instead, 

^ Nsi/2D~1 

,sl,i 1 9rj,sl \9fj,sl) 1 9fj,sl,i / ^ l^g,! ^Vj ,sl,i,g i 



9=0 

[Nsi/2D-1] 

E 1^-1' = ^' (37) 

5=0 

where D = 1,2 and the operator gf-^si,i acts onto the spin and si effective lattice 
associated with the A^j\ = 1 configuration state. For each sub-configuration state the 
real-space coordinates fj of the si bond particles are the same as for the A^^^^ = 
configuration state except for the corresponding d, I line. For that line they are shifted 
by 2fj^i as given in Eq. (I36l) . In the above equation we denote by d" and /" the summation 
indices in the expression for the operator aff.^si,i,g given in Eq. ( I37|) to distinguish them 
from the constant d, I line indices d and /. 

The two-site one-bond operators bi^^si,i,d",i",g and b],siid"i" g involved in Eq. (1371) 
have for si bond particles whose real-space coordinates do not belong to d, I, I' lines the 



9rj,sl = 


V2BY*'"" 


0'fj,sl,i,g 


d" = l /"=±1 
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same expression for all 2D = 2,4 sub-configuration states and read, 

'^r,sl,i,d",l",g ^ '^f,sl,d",l",g ' (38) 

where bt^^ ^„ ;„ is given in Eq. f l27|) . In turn, when fj belongs to the dj,0 = dj line 
such operators are given by, 

b^ 



(-1) 



d+d" 



V2 



2 + ■5^_^9 



d",l" 



^^^I'lu'^n 



2 ^ ^^^^l".v'^^- 
~ ^d"4 2?^d/ , 



d"j" 



(39) 



where we recall that 1 = 2 and 2 = 1. Finally, for the si bond particles whose real- 
space coordinate r,- belongs to the d,l,l' lines such that /' = ±1 the two-site one-bond 
operators read. 



bl 



f,sl,i,d",l",g 



V2 ^ 


'1 
2^' 


z 

' 'd",l" 


'1 

2 + '^^+^^d 


II ill^^i 


^ ^d",l 



d",l" 



^d",d^l"-l''^'^d,l 



(40) 



4-5. More about A^j\ = 0, 1,2 configuration states 

The si fermion momentum occupancies of Ref. [2] of excited states generated by 
application onto x > and m = ground states of two-electron operators involving 
(i) creation or annihilation of two electrons with the same spin projection and (ii) 
spin-triplet and spin-singlet excitations are described by a superposition of the set of 
[A^si -|- 2] [Nsi + l]/2 configuration states with A^^^^ = 2 unoccupied sites studied in the 
following. Indeed, there are [A^^i + 2] [A^^i + l]/2 different A^^^ = 2 configuration states 
corresponding to the possible [A^'^i -|- 2] [A^^i + l]/2 positions of the two unoccupied sites 
in the si effective lattice. According to Eq. ([9]), here A^^i + 2 = [N^J2 + Ss] where 
Ss = I for Ns2 = and Sg = for Ns2 = 1 is the number of sites N^^_^ of that lattice 
whose number of unoccupied sites reads A^j\ = [25*5 -|- 2Ns2] = 2. For such states the 
number N^^ = [2Nsi +A — 2Ss] of sites of the spin effective lattice is even. There are two 
sets of [Nsi + 2] [Nsi + l]/2 configuration states which contribute to different N^^ = 2 
energy eigenstates. Those are the A^j\ = 2, Ss = 1, Ns2 = spin-triplet configuration 
states and A^j\ = 2, 5"^ = 0, Ns2 = 1 spin-singlet configuration states, respectively. 

Since it is assumed that for the square lattice [2iVsi]^'^ is an integer number, the 
notations N^^_^ and A^^^ do not mean that the numbers iVa^-^ and A''^^ are integers. 
^asi ~ ^si + 2 and A"^^ = [2Nsi + 4 — 2Ss] are integer numbers. In turn, the integer 
numbers closest to A"a^^ and Na^ are the mean value of the number of sites of each 
row and column of the si and spin effective lattice, respectively. Our study takes into 
account the exact number of sites of such rows and columns. 
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Concerning the above two types of A^^^^ = 2 configuration states, the two unoccupied 
sites can either correspond to two independent +1/2 spinous {Sg = 1 and Ns2 = 0) or two 
spinous out of the four which are part of the s2 fermion {Sg = and Ns2 = 1). We recall 
that alike an independent spinon, the spin-neutral four-spinon occupancy configuration 
of such a s2 fermion remains invariant under the electron - rotated-electron unitary 
transformation [2]. 

For the N^ — )■ cxd limit considered here for the model on the square lattice the 
correct physics is achieved if one considers that out of the [2A^si + 4 — 2Ss] sites of the 
spin effective lattice, 2Nsi sites correspond to a square lattice and the position of the 
two extra sites which contribute to the si effective lattice is defined as for the N^^ = 1 
configuration states. For the model both on the ID and square lattice and Ss = and 
thus Ns2 = 1 the number of sites of the spin effective lattice is given by A^^ = 2Nsi + 4. 
However, two of the four extra sites do not contribute to the si effective lattice. The 
point is that for A^^^^ = 2 configuration states with Sg = and Ns2 = 1 when the si bond 
particles moves around in that lattice, they use as unoccupied sites only two out of the 
four sites of the s2 fermion. Two of these four sites of the spin effective lattice do not 
belong to the si effective lattice, alike the sites of the ?7-spin effective lattice associated 
with the original rotated-electron unoccupied sites do not belong to the spin effective 
lattice. 

The spin SU{2) symmetry imposes that the si effective lattices of the N2 = 2 
configuration states such that (i) Sg = I for Ns2 = and (ii) Sg = for Ns2 = 1, 
respectively, must be identical. In either case each of the two extra sites obeys the 
boundary conditions given in Eq. ( 1M|l which for the present A^^^ — )■ 00 limit yield the 
correct physics. For the square lattice the row and column periodic boundary conditions 
are similar to those of the N^-^ = 1 configuration states. For the ID lattice periodic 
conditions are also used. 

Symmetry imposes that for gf-primary partitions of A^j\ = 2 configuration states the 
two unoccupied sites are associated with kink-like and anti-kink-like link configurations, 
respectively. Moreover, symmetry implies that only sub-configuration states for whose 
5f-primary partitions both the kink and anti-kink link configurations associated with 
the two unoccupied sites, respectively, correspond to rows or columns are allowed. 
It follows that alike for A^j\ = 1 configuration states, a N^-^^ = 2 configuration 
state involves 2D = 2,4 sub-configuration states. In the present case there are 
two horizontal sub-configuration states for whose f^-primary partitions the unoccupied 
sites of real-space coordinates fj^ and fji are associated with row kink-like and anti- 
kink-like link configurations and row anti-kink-like and kink-like link configurations, 
respectively. Furthermore, there are two vertical sub-configuration states for whose g- 
primary partitions the unoccupied sites with such real-space coordinates are associated 
with column kink-like and anti-kink-like link configurations and column anti-kink-like 
and kink-like link configurations, respectively. 

It is useful to classify the A^^^ = 2 configuration states into three main classes: 
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1) A^^^ = 2 configuration states whose two extra sites do not belong to the same 
row and column. Such states occur for the model on the square lattice but do not exist 
for the ID lattice. For the corresponding (yf-primary partitions the occupied sites are the 
extra sites, which correspond to kink-like or anti- kink-like link configurations of the same 
form as those of N^^ = 1 states. Indeed, each unoccupied site behaves independently and 
can move around in the spin effective lattice by means of the elementary steps discussed 
above for the A^^^^ = 1 configuration states. For all gf-primary partitions one unoccupied 
site is associated with a kink-like link configuration and the other with an anti-kink-like 
link configuration, respectively. The si bond-particle operators are defined as for the 
Ngi = 1 configuration states by Eqs. ( l371) - (H0l) but with parallel d,+l and d,—l lines 
and the index i corresponding to the four sub- configuration states defined above for the 
iV^'\ = 2 configuration states. 

2) N^i = 2 configuration states whose two extra sites belong to the same row or 
column but are not nearest neighboring sites. For the ID lattice only states of this 
type and of type 3 exist. Let the two extra sites have real-space coordinates fj^ and 
fj' . Then for the ^f-primary partition of one of the two horizontal (and vertical) sub- 
configuration states the link occupancy involving sites belonging to the row part (and 
column part) located between the sites of real-space coordinates Vj^ and fj'^ has opposite 
sign relative to those of the surrounding rows (and columns) whereas the link occupancy 
involving sites belonging to the row part (and column part) located between the sites 
of real-space coordinates fj' and fj^ has the same sign as those of the surrounding 
rows. In contrast, for the (jf-primary partition of the other horizontal (and vertical) 
sub-configuration states the sign of the link occupancies on the two above row parts 
(and column parts) is the opposite. For (yf-primary partitions the two unoccupied sites 
of real-space coordinates fj^ and fj' behave as kink or anti-kink mobile domain walls 
whose link structure is shown in Figs. 7 and 9. The si bond-particle operators are 
defined as for the N!!-^ = 1 configuration states by Eqs. fl37|) - P0|) but with the operators 
whose real-space coordinates correspond to link occupancy configurations with opposite 
sign relative to those of the surrounding rows or columns referring to the d, j line parts 
defined above instead of to the whole d, j line. 

3) ^si — 2 configuration states whose two extra sites belong to the same row 
and column and are nearest-neighboring sites in the spin effective lattice. In that case 
there is both a horizontal and a vertical sub-configuration state whose extra-site row and 
column link occupancies have opposite sign relative to those of the surrounding rows and 
columns, respectively. In turn, for the other horizontal and vertical sub-configuration 
states, respectively, all link occupancies have the same sign and the two-unoccupied-site 
link occupancy configuration of the corresponding gf-primary partitions are identical. A 
small sub-domain of such a link configuration including the two-unoccupied sites of the 
spin effective lattice is shown on the right-hand side of Fig. 4. These configuration states 
can be generated from those of type 2. This is achieved by moving both unoccupied 
sites through their row or column, one against the other, until they become nearest- 
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neighbor sites. For the square (and ID) lattice, the g-phmarj partitions of two (and 
one) out of the four (and two) sub- configuration states such a motion and resulting 
" collision" leads to the full cancellation of the corresponding kink-like and anti-kink-like 
link occupancy configurations. Such a unoccupied-site kink and anti-kink annihilation 
has similarities with a particle - anti-particle annihilation process. The corresponding 
"collisions" occur more often for a unoccupied site moving around in the spin effective 
lattice of the ID lattice. Indeed, then there is a single chain. In contrast, for the model 
on the square lattice there are both horizontal and vertical elementary steps so that the 
two unoccupied sites can move independently around in the spin effective lattice without 
colliding provided that their motion involves configuration states of type 1, which do 
not exist in ID. For the two above sub-configuration states of the model on the square 
lattice whose extra-site row or column link occupancies have opposite sign relative to 
those of the surrounding rows or columns, respectively, the si bond-particle operators 
are defined as for the N^^ = 1 configuration states by Eqs. (l371) - (H0l) . For the other two 
sub-configuration states whose extra-site row or column link occupancies have the same 
sign as those of the surrounding rows or columns, respectively, the si bond-particle 
operators are defined alike for the N^^ = configuration states by Eq. ( l24l) . 

The number A^j\ of unoccupied sites is a good quantum number for the model on the 
square lattice in the one- and two-electron subspace so that the si fermion momentum 
occupancies of Ref. [2] of the states that span such a subspace are described by a suitable 
superposition of configuration states with constant number of unoccupied sites in the 
si effective lattice. For instance, the si fermion momentum occupancy of a N!!-^ = 2 
spin-triplet (and spin-singlet) excited state is described by a suitable superposition of 
the [N,i + 2] [A^,i + l]/2 A^,'^ = 2 configuration states with S, = 1 and N,2 = (and 
with Ss = and Ns2 = 1-) Application of a si bond-particle creation operator of real- 
space coordinate fj onto such a ground state gives zero for configuration states of the 
above types 1 and 2: only the two sub-configuration states of a configuration state of 
type 3 whose link configurations have all the same sign and the real-space coordinate of 
the central unoccupied site on the right-hand side of Fig. 4 coincides with that of the 
applied operator are transformed onto the N^^ = configuration state with one more 
si bond particle than the initial state. 

Moreover, application of a si bond-particle creation operator onto N^^ = 0, 1 
configuration states gives always zero. In turn, application of a si bond-particle 
annihilation operator whose real-space coordinate is that of an occupied site of the 
si effective lattice of a N^^ = 0,1,2 configuration state transforms it into a state with 
two more unoccupied sites than the initial state. 

Finally, in Appendix A it is confirmed that within the description used in the 
studies of Ref. [2] one has that for U/4t > and the subspaces where the si bond- 
particle operators of the one- and two-electron and N^^ = 0,1,2 subspace are defined, 
upon acting onto the si effective lattice such operators anticommute on the same site 
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of that lattice, 

{9l^,sv9rj,sl} = 1 ; (4,^1' 4:^1^ = {9rj,au, 9r„sl} = , (41) 

and commute on different sites, 

[9t^,sv9ry,sl] = [4,sl' ^'-^l^ " l9r„sl, 9ry ,l] = , (42) 

where j ^ j'. That algebra confirms that the si bond-particle operators are hard-core 
like. 

5. Concluding remarks 

For the Hubbard model on the square lattice in the one- and two-electron subspace the 
number NJ!^ of unoccupied sites of the si effective lattice is a good quantum number. 
Such a subspace is spanned by the states with N^^ = 0, 1, 2, whose spin configurations 
are investigated in this paper. In the studies of Sections 3 and 4 a change of gauge 
structure |22j is considered so that the real-space coordinates of the sites of the si 
effective lattice correspond to one of the two sub-lattices of the spin effective lattice. 
That change is fulfilled for the A'^^^ = configuration state for which the gauge structure 
occurs. That leads to two alternative definitions of the si effective lattice for both that 
state and the related A^^^^ = 1,2 configuration states: As discussed in this paper, the 
occupancy configurations of the latter states can be generated from those of the N!!^ = 
configuration state with the same number of si bond particles. 

The N^^ = configuration states studied here generate the spin degrees of freedom 
of the X > and m = ground states and their two-electron charge excited states of 
the Hubbard model on the square lattice in the one- and two-electron subspace. That 
quantum problem refers to the square- lattice quantum liquid of c and si fermions of 
Ref . [2] . In turn, the spin configurations of the excited states generated from application 
onto these ground states of one-electron operators and two-electron operators other 
than the charge operator are suitable superpositions of the A^^^ = 1 and N^^ = 
or N!}-^^ = 2 configuration states, respectively, studied in this paper. For the A^^^^^ = 1,2 
states the unoccupied sites move around in the spin effective lattice through well-defined 
elementary steps. For A^^'*^ = 2 configuration states the two unoccupied sites may move 
through the same row or column, one against the other, until they become nearest- 
neighbor sites. For the model on the square (and ID) lattice, for two out of four (and 
one out of two) (/-primary partitions defined in this paper such a motion and resulting 
"collision" leads to the full cancellation of the corresponding kink-like and anti- kink- 
like link occupancy configurations, similarly to a particle - anti-particle annihilation 
process. Such collisions occur more often for an unoccupied site moving around in the 
spin effective lattice of the model on the ID lattice, which corresponds to a single chain. 
As confirmed in Appendix A, the spin configurations obtained from such a cancellation 
play an important role in the sl bond-particle operator algebra. 

Moreover, in this paper suitable sl bond-particle operators have been constructed, 
which upon acting onto the sl effective lattice obey a hard-core like algebra. Such a 
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property is valid for the subspaces in which the si bond-particle operators act onto. It 
plays an important role in the related studies of Ref. [2] . There corresponding si fermion 
operators are generated by means of an extended Jordan- Wigner transformation from 
the si bond-particle operators introduced in this paper. Concerning previous studies on 
the large-f/ Hubbard model and t — J model on a square lattice involving for instance 
the slave particle formalism p!6| [T8l [22] or Jordan- Wigner transformations [12], the 
crucial requirement is to impose the single occupancy constraint. Here that constraint 
is naturally implemented for all values of U/4t > 0, since the spins associated with the 
spin SU{2) state representations refer to the rotated electrons of the singly occupied 
sites. Moreover, for the above schemes the spinless fermions arise from individual spin- 
1/2 spins or spinous. In contrast, within the extended Jordan- Wigner transformation 
performed in Ref. [2] the si fermions emerge from the spin- neutral two-spinon composite 
si bond particles studied in this paper. 

The studies of Ref. [2] confirm that the results obtained here concerning the 
Hubbard model on the square lattice in the one- and two-electron subspace are useful for 
the further understanding of the role of electronic correlations in the spin- wave spectrum 
observed in the parent compound La2Cu04 [Ij. A system of weakly coupled planes, each 
described by the square-lattice quantum liquid of Ref. [2], perturbed by the effects of 
intrinsic disorder is expected to be the simplest realistic toy model for the description of 
the role of correlations effects in the unusual properties of the cuprate hight-temperature 
superconductors [HI [151 UHl E] • 

I thank Nuno M. R. Peres for discussions and support in the figures production. 
I also thank Miguel A. N. Araiijo, Daniel Arovas, Pedro D. Sacramento, and Maria J. 
Sampaio for discussions and the support of the ESF Science Program INSTANS and 
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Appendix A. Hard-core character of the si bond-particle operators 

The goal of this Appendix is to confirm the validity of the relations provided in Eqs. (|4T!) 
and (H2l) . In order to probe such relations, we consider without any loss of generality that 
the initial state is a N^^ = 4 configuration state. Indeed and as discussed in Subsection 
4-5, application of a si bond-particle creation operator onto N^-^ = 0, 1 configuration 
states gives zero. Application of two si bond-particle creation operators onto N^-^ = 2 
configuration states gives zero as well. The studies on NJ!^ = 2 configuration states of 
that subsection can be straightforwardly generalized to N^^ = 4 configuration states. 
That assures that application onto those of four-site two-bond operator may not give 
zero. Fortunately, the final results reached here apply to any configuration state with a 
finite number Nj}^ of unoccupied sites in the si effective lattice, yet are simpler to derive 
for the N^^ = 4 configuration states. 

Sub-configuration states of type 3 considered in Subsection 4-5 for N^-^ = 2 also 
exist for N^^ = 4. Within those with A^^^^ = 4 unoccupied sites, the importance of 
sub-configuration states whose link configurations have all the same sign stems from 
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those being the only ones for which apphcation of si bond-particle creation operators 
does not give zero. For other A^j\-finite sub-configuration states the si bond-particle 
operators have slightly more involved expressions given in Section 4. Fortunately, the 
use of the operator expressions suitable for sub-configuration states of type 3 leads to the 
same final results. Therefore, for simplicity we consider here the latter si bond-particle 
operators whose expressions are the same as for the N^-^ = configuration state studied 
in Section 3. 

All the results of this Appendix apply to the model on both the square and 
ID lattices. The ID expressions are readily obtained if one considers in the general 
expressions given below only d = 1 contributions and terms, together with the choice 
D = 1 in the Z)-dependent quantities. In order to reach our goal, let us recall that the 
rotated-electron operators are related to the original electron operators as follows, 

cl. = V^cl.V^ (A.l) 

where V is the electron - rotated-electron unitary operator. Straightforward 
manipulations based on Eqs. ([5])- ([7]) then lead to the following algebra for the related c 
fermion operators [2], 

{ffj,c , /,=*.,c} = ^j,j' ; ^fr,,c ' fr,,,c} = {frj,c , fr.,,c] = , (A.2) 

c fermion operators and rotated-quasi-spin operators, 

[4c'4] = [^.>-4] = 0' (A-3) 



J 





p" 


{^'^r^^ 


= 1, Kv^?> = o 


[4,^%]- 


= <5,,,.2g;?S [g±,g± 



and rotated-quasi-spin operators, 

p=l,2,3, (A.4) 

(A.5) 
1 = 0. (A.6) 

J 

Hence the rotated-quasi-spin operators q^, anticommute on the same site and commute 
on different sites. 

Moreover, combining the expressions given in Eq. (|6]) with the algebraic relations 
provided in Eqs. ( 1A.4l) -( lAl6l) one readily finds that, 

{s+,sz} = l, {s±,s±} = 0, (A.7) 

Kv^.-J = Kv4] = 0' (A.8) 

for i ^ i' and, 

[4'4]=±^«'^.V (A-9) 

It follows that the spinon operators s^ anticommute on the same site and commute 
on different sites. Consistently with the rotated-electron singly-occupied site projector 
Uf^-^c appearing in the expression of the spinon operators s^ and s^ provided in Eq. 
(jH]), their real-space coordinates fj can in the present N^ — )■ oo limit be identified with 
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those of the spin effective lattice. Therefore, the corresponding operator index values 
j = 1, ..., Na^ are chosen to be those of that lattice. 

For the sub-configuration states of type 3 the operators gt ^^ (and gfj,si) which 
create (and annihilate) a si bond particle at a site of the spin effective lattice of real- 
space coordinate fj have the general form given in Eq. (|2^ both for the model on the 
ID and square lattices, where the absolute value \hg\ of the coefficients hg decreases 
for increasing link length C,g and obeys the normalization sum-rule fl25l) . Hence the 
expression of the si bond-particle operator gt ^-^ involves the operators at _,^ and arj,si,g 
of Eq. i^^, which create and annihilate, respectively, a superposition of 2D = 2, 4 two- 
site bonds of the same type and bl^^^, and bf^si,d,i,g are two-site one-bond operators 
whose expression is given in Eq. fl27j) . 

In order to confirm the vahdity of Eqs. (jH]) and (|42l) let us use Eqs. (l24l) -( !27l) 
to rewrite the anti-commutation relations of Eq. ( HTi) in terms of anti-commutators of 
two-site one-bond operators as follows, 

(4,^1' 3rl,sl} = J2J2K f'a' H^+rO^,sl Al,9^ K^.+rO ^„sl,d' ,1' ,g'} , (A.IO) 

d,l,g d',l',g' 

d,l,g d',l',g' 

{9rj,sl,grj,sl} = 2_^ 2_^ ^9 V{^r^+rO,,sl,(i,Z,g)^r^+rO j,,sl,d',r,<;'} ; (A.12) 

d,l,g d',l',g' 

where for simplicity we used the abbreviated summation notation, 

D N,i/2D-1 

E-EE E ■ (A.13) 

d,l,g d=l l=±l 9=0 

Moreover, by the use of the same equations the commutation relations of Eq. (142|) can 
be expressed in terms commutators of two-site one-bond operators. That leads to, 

[9l^,sv9ry,sl] =J2J2K V [^J- +.-0,,sl,d,L<,' ^,+rO_,„sl,d',Fy] , (A.14) 

d,l,g d',l',g' 

[9r^^,sv9^^,,sl]= Z^ Z^ KK'\-^rj+rO^,sl,d,l,g'^rj/+r^, ^nsl4',l',g'h (A.15) 

d,l,g d',V,g' 

[9rj,sl,9r^,,sl\=2_^ Z_^ ^9^9'\Prj+rli,sl,d,l,g-,^r^,+r'^,^,,sl,d',V,g']^ (A. 16) 

d,l,g d',l',g' 

where j ^ j'. 

According to the studies of Subsection 3-3, four rules follow from the definition of 
the subspace where the operators of Eqs. (I24l)-(l271) act onto. The evaluation of the 
anti-commutators and commutators of the two-site one-bond operators on the right- 
hand side of Eqs. (1A.10l) -( p^.l2l) and ( lA.14l) - (1A.16ll . respectively, relies on both such 
rules and the algebra given in Eqs. ( 1A.7I) - (1A.9I) of the spinon operators s^ and sp, of 
Eq. dnj, which are the building blocks of the two-site one-bond operators of Eq. (|5j 
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Fortunately, according to Eqs. (IA.7p -( 1Xl9|) the spinon operators s^ obey the usual 
algebra: They anticommute on the same site of the spin effective lattice and commute 
on different sites. 

The two-site one-bond operators of Eq. (!27|) can be rewritten as, 






:-i 



,d-i 



V2 



l + s; 



'r,- 



\^'^ 



(A.17) 



and h 



ri,r2 



bl ^ ) where recalling that the real-space coordinate of their link centre 





reads r = rj + r^^ the real-space coordinates ri and f^2 are given by. 



ri 



^7 + r,j - rli ; f2 = r,' + r,, ,, + f|, ,, . 



' 3 ' ' d,i 'd,n '^~'j' I 'd',i' ' 'd',i'- (A. 18) 

The evaluation of the anti-commutators and commutators of the two-site one-bond 
operators on the right-hand side of Eqs. ( 1A.10P -( 1A.12I) and ( ]A.14[) -( lA.16l) . respectively, 
then relies on straightforward manipulations based on Eqs. ( ]A.7P -( lAl9l) and on the four 
rules given in Subsection 3-3, which define the subspace that the operator algebra under 
consideration refers to. For the two general anti-commutators needed to evaluate the 
two-site one-bond operators on the right-hand side of Eq. (lA.lOp we find the following 
expressions. 



{bk 



^2? '^ri,r2 



} 



E(ia 



«=1,2 



S^ S^ 
Ti r-. 



2 ' 



- + si 
2 + ^-1 



1 



'- + s^ 



'- + s^ 



(A.19) 



{bi 



■,} 



4 = 1,2 ^ 



1 

2+4 



s^ s^ i — 



o + 4- 

2 ' 



} 



(A.20) 



where 1 = 2, 2 = 1, ri ^ fr , r2' , and f2 ^ rv,r2'- Indeed, according to the second rule 
for the subspace where the two-site one-bond operators act onto only general operators 
6t ^2&^^,,r2/ which do not join sites or join both sites of the spin effective lattice lead to 
wanted and physical spin configurations. 

Moreover, in the initial configuration that the four-site two-bond operators of the 
relation ( 1A.19|) act onto one has according to the first and third rules that the two sites 
of real-space coordinates ri and r2 either (i) are occupied by two independent -1-1/2 
spinous or (ii) are linked by a bond. In turn, the anti-commutator of Eq. (]A.20|) is of 
the form of those on the right-hand side of Eq. ( JA.IOI) so that in Eq. (IA.18I) one has 
that j = j' yet f^i ^ r*J y and the restrictions imposed by the forth rule must be taken 
into account. Since we find that 6t ^b^^i^r^, = bf^,^f^,bt^ ,^^ and each of such operators is 
given by one half the operator on the right-hand side of Eq. (JA.201) . when the two sites 
of real-space coordinates fi and r2 (and ri/ and r*2') are linked by a bond (and occupied 
by two independent -1-1/2 spinous) one must consider both initial configurations where 
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the sites of real-space coordinates fi/ and 7^2' (and fi and r2) are (i) linked by a bond and 
(ii) occupied by two independent +1/2 spinons. In turn, when the two sites of real-space 
coordinates fi and r2 (and fii and r2') are occupied by two independent +1/2 spinons 
(and linked by a bond) one must consider only initial configurations where the sites of 
real-space coordinates ri' and r2' (and n and r2) are linked by a bond (and occupied 
by two independent +1/2 spinons). 

It then follows from analysis of the operator expression on the right-hand side of 
Eq. ( ]A.19|) that when in the initial configuration the two sites fi and r2 are occupied 
by two independent +1/2 spinons the operator term X]j=i 2[V2](l/2 + s^Y transforms 
that configuration onto itself whereas the remaining operator terms give zero. In turn, 
when in the initial configuration the two sites fi and f^2 are linked and correspond to an 
one-bond configuration the operator terms J2i=i 2[l/2]{(l/2 + spY — stsZ(l/2 + s^J} 
transform that configuration onto itself whereas the remaining operator terms give zero. 
On the other hand, when acting on the above initial configurations the operator on the 
right-hand side of Eq. (1A.20I) gives always zero so that when acting onto the subspace 
that the operators of Eqs. fl24|) - fl27j) refer to the anti-commutators provided in Eqs. 
( JA.191) and ( JA.201) simplify and are given by, 

H,,r,M,r2} = 1; {bl„,,,h'A'} = 0^ (A-21) 

where fi ^ fi/ , f2' and r2 ^ fi/ , f2' . 

Next concerning the two general anti-commutators needed to evaluate the two- 
site one-bond operators on the right-hand side of Eq. flA.llI) we find the following 
expressions, 

{bkrAur.} = -2 (^ + 4) [I + 4) 4^.^.-E {I + 4.) 4V. , (A.22) 



Sr.Sr., 



where as above 1 = 2, 2 = 1, ri 7^ n', "^2', and r2 7^ ^i', "'^2'- 

It follows from analysis of the operator on the right-hand side of Eq. ( IA.22|) that 



according to the first and third rules when in the initial spin configuration the sites of 
real-space coordinates fi and r2 are both (i) occupied by independent +1/2 spinons and 
(ii) linked by a bond, application of that operator onto such a configuration gives zero. 
In turn the anti-commutator of Eq. ( JA.23P is of the form of those on the right-hand side 
of Eq. (lA.lip so that in Eq. (]A.18|) one has that j = j' yet r^^ 7^ r|, j, and then the forth 



rule applies. Since we find that bl -.bl ^ = bl -, bt -, and each of such operators is 
given by one half the operator on the right-hand side of Eq. (JA.23P , concerning the latter 
operator when the sites of real-space coordinates fi and r*2 (and fi/ and r2') are linked 
by a bond one must consider both initial configurations where the sites of real-space 
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coordinates fii and r2/ (and r*i and 7^2) are (i) linked by a bond and (ii) occupied by two 
independent +1/2 spinons. In turn, when the two sites of real-space coordinates fi and 
r2 (and fir and f^2') are occupied by two independent +1/2 spinons one must consider 
only initial configurations where the sites of real-space coordinates fir and r2/ (and fi 
and r2) are linked by a bond [2]. It follows then from analysis of the operator on the 
right-hand side of Eq. flA.23|) that application of it onto any of such spin configurations 



gives zero. 

A similar analysis for the two general anti-commutators needed to evaluate the 
two-site one-bond operators on the right-hand side of Eq. (1A.12|) leads to, 



{b,,,,,, b,,,,,} = -244 Q + si}j Q + si}j - Y^ 44 Q + st^ , (A.24) 



^f^^r / 



where as above 1 = 2, 2 = 1, ri 7^ ri/, r2/, and r2 7^ f^, r2'. 

Again analysis of the operator on the right-hand side of Eq. ( 1A.24I) reveals that 
according to the first and third rules when in the initial spin configuration the sites of 
real-space coordinates n and r2 are both (i) occupied by independent +1/2 spinons and 
(ii) linked by a bond, application of that operator onto such a configuration gives zero. 
On the other hand, the anti-commutator of Eq. (IA.25P is of the form of those on the 



right-hand side of Eq. (JA.lip so that in Eq. ( lA.lSp one has that j = j' yet f^i 7^ f^, ;, 
and then the forth rule applies. Since we find that bTf^fJ)^^^,^^, = &r^,,r2/^ri,r2 and each of 
such operators is given by one half the operator on the right-hand side of Eq. (1A.23I) . 
concerning the latter operator when the sites of real-space coordinates fi and r2 (and f^ 
and r2/) are occupied by two independent +1/2 spinons one must consider both initial 
configurations where the sites of real-space coordinates fr and r2' (and fi and r2) are 
(i) linked by a bond and (ii) occupied by two independent +1/2 spinons. In turn, when 
the two sites of real-space coordinates fi and f2 (and fy and r2/) are linked by a bond 
one must consider only initial configurations where the sites of real-space coordinates 
fir and r2' (and fi and r2) are occupied by two independent +1/2 spinons [2]. Analysis 
of the operator on the right-hand side of Eq. f lA.23|) then reveals that application of it 
onto any of such spin configurations gives zero. 

It then follows from the above results that when acting onto the subspace that the 
operators of Eqs. f l2^ -f l27|) refer to the anti-commutators provided in Eqs. flA.22|) and 
( 1AI23|) and Eqs. ([Qill and ( lAl25l) simplify and read. 



(4,,^,, 4,,-,} = (4,,^,, bl^,,J = ; {b,,,,,, 6^,,-,} = {6,-.,,^,, b,^„,J = , (A.26) 
where r 1 7^ r 1/ , r2/ and r2 7^ r 1/ , r2' • 
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The use in Eq. (lA.lOp of the anti-commutators of Eq. (lA.2ip with the two-site 
one-bond operators related to those of Eq. fl271) by the expressions provided in Eqs. 
(IXTtI) and (lAlSJl leads to, 

{9lsv9r.,si} = Yl \^f = 2DJ2M = 1, (A.27) 

d,l,g g 

which is the first relation of Eq. ( 14T1) . To perform the summation of Eq. ( 1A.27I) the 



sum-rule (l25l) was used. Furthermore, the use of the anti-commutators of Eq. (JA.26P in 
Eqs. (lA.llI) and ( 1A.12I) leads to the remaining relations of Eq. (1411) . 



The evaluation of the commutators of Eq. fl42|) by the use of the expressions given 
in Eqs. flA.14|) - flA.16|) is much simpler. First it is simplified by the property that two- 



site one-link bonds belonging to si bond-particle operators with different real-space 
coordinates are always different. Second the evaluation of such commutators also relies 
on straightforward manipulations based on Eqs. (IA.7I) - (1A.9I) . which lead directly to, 

[^kr,^ K.r.>] = [bkr.^ ^^^rj = [^n/., ^„rj = , (A.28) 

for r 1 7^ r 1/ , r2/ and r2 7^ r 1/ , r2' . 

Finally, the use in the expressions of the si bond-particle operators of Eq. ( 142|) of 

the commutators of Eq. (1A.28P with the two-site one-bond operators related to those of 

Eq. ((27]) by the expressions provided in Eqs. f lA.171) and f lA.lSD in Eqs. f lA.14D - flA.16D 

leads readily to the commutation relations provided in Eq. fH2l) . 
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